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THE FLOW OF HEAT THROUGH PLATES! 
By R. RueEpy? 


Abstract 


The speed is calculated with which the steady flow of heat is established in a 
slab of uniform temperature after one boundary plane has been suddenly brought 
to a higher tem ture, or when the temperature of both planes is changed. In 
both cases the flow of heat may be expressed by means of simple theta functions, 
and the law of approach to the steady state may be used for determining the 
diffusivity of the material. When one boundary plane undergoes sinusoidal 
variations in temperature while the other is maintained at a constant level, a 
finite thickness is found for which, in the steady state, the heat flowing in or out 
during one half-cycle reaches its highest value. 


The Steady State in the Plate Method 


When the “plate’’ method is used for measuring the heat conductivity of 
solid bodies, one side of a slab is more or less suddenly brought into contact 
with a source of heat at constant temperature, while the other side is kept 
cool. Heat begins to flow through the slab, and when a steady state is reached, 
the coefficient of heat conductivity, k, is determined. Losses which would be 
caused by radiation, conduction and convection through the rest of the 
boundary surface are prevented with the aid of a guard ring surrounding the 
slab. How closely the heat conductivities may be ascertained depends on the 
accuracy with which the thickness of the plate, the temperatures of the two 
walls, and the heat energy, Q, furnished to one side or drawn from the colder 
end, may be measured. The coefficient k is determined from the formula 

Q(watts) Xb(cm.) 
4.184A (sq. cm.) X (0: —0%) deg. C. 


where v, — vis the temperature difference maintained between the two bcundary 
planes. As it is not practicable, in general, to measure the thickness to within 
less than 0.1 mm. it is desirable to test slabs at least 5, or better, 10 mm. thick 
in order to keep the possible error below 2%. For a temperature difference 
of 15° C., an interval over which k may be considered as being constant, an 


1 Manuscript received May 12, 1932. 


Contribution from the National Research Laboratories, Ottawa, Canada. This paper was 
read at the meeting of the Royal Society of Canada, May 1932. 


2 Junior Research Investigator, National Research Council, Ottawa. 


k(c.g.s.) = 





578 CANADIAN JOURNAL OF RESEARCH 


error of 0.2° C. means an uncertainty of 1.3%. The energy, Q, furnished, 
computed from the current measured to within 0.002 amp. and the voltage 
read to within 0.04 volt, is known with an accuracy of between } and 1%, so 
that the uncertainties entering into the value of & will not be below +2% (2). 
The precision may be improved by using thicker plates, but when many 
measurements have to be made the question arises as to how much longer it 
will take for the steady state to be established. While it is also possible to 
increase the temperature difference between the hot and the cold plane, the 
changé of k with temperature may introduce new errors. 

When at the beginning of the experiment the temperature in a slab of thick- 
ness b, specific heat c and density s, is a function, f(z), of the distance z from 
the cold wall of the slab, and is equal to zero at time ¢ equal to zero, while the 
ends are kept at the temperatures Fo(/) and F,(/) during the interval from 
t= to i=1, then the temperature v at any distance z is at any instant ¢ repre- 
sented by the following formula (1, p. 69). 

oo * ‘ ' a n 
ee Fa [serine Fades * f me Ra) —(—1) roa | 
o o 


nel 
where x= = is the diffusivity of the material. 


In the simplest case the temperature is constant throughout the thickness, 
and one side of the slab is maintained at this temperature which may, therefore, 
be chosen as the zero point so that f(z) becomes equal to zero. The first 
term in the brackets in the general equation vanishes, so that F,(z) = V being 
the constant temperature at the hot side, the relation reduces as follows: 
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Integrating and taking into account the fact that 

et _ 2 +46 3x 

sin b z—4 sin b 3% sin b 3 
the equation becomes 

2 = (-1)" «= 
ye v ($425 € sine Fe) 
b = n=1 n 

The heat C flowing through unit area at a distance z and time ¢ is obtained by 
differentiating v with respect to z and multiplying by &. 
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n=l 


C=-—k ~ ote, q) 


-# 
where @ is the symbol for the theta function and g=e ® _ The heat 
directly measured and flowing into the slab from the hot side z=) is therefore 
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given by the expression 


V 
Cy=—k ~ @3(0, g) = C%@,(0, 9), 


whereas the heat penetrating into unit area of the cold wall is 


Gah onto q) =C°@x(0, 9), 


where C° is the heat flow ( - rr in the steady state. The values of these 


theta functions are listed in convenient recently issued tables (3). As an 
illustration of the rate at which the steady state is established Table I shows 
the times in seconds in which certain fractions 0; and @ of the final heat flow 
are reached in the case of a slab of sulphur 3.14 cm. thick (conductivity , 0.0005; 
specific heat, 0.17; density, 2.0; and therefore x, 0.00147), or slightly com- 
pressed spruce of 16% moisture content, the flow being measured parallel to 
the grain (specific heat, 0.33; density, 0.4). Under the conditions stated the 
coefficient of heat conductivity is the same for sulphur and spruce but on 
account of the smaller value of cs a higher value for x (0.0038) results for spruce, 


TABLE I 
TRANSIENT STATE OF HEAT FLOW 


Time in sec. 


xt exp. (—xF+) 63 b=x cm. 
Sulphur Spruce 


1041 404 

1327 515 

1585 

1843 

2184 

2578 1000 

2905 1127 

3476 1348 

| : 3660 1420 
4633. | 1797 

| 5177 2008 | 


NOM Ue WWD De 
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and less time is needed for reaching the steady state in which 0;=@)=1. 
The times given for sulphur also apply to wood fibre (k, 0.00015; c, 0.31; 
s, 0.33; x, 0.00146). For any other material they may be found by dividing 
the figures of. the first column by x. As is to be expected, the rate at which 
heat flows at the beginning from the hot side into the slab is larger than the 
amount flowing in the steady state because the material must be heated up. 
Near the cold side, on the contrary, there is practically no heat flowing until 


wal or xt, in the present case, exceeds a value of about 0.23; after it has become 


equal to 0.3, the flow is 0.002 of the final value Co, at 0.6 equal to 0.074 C°, 
at 1.0 equal to 0.3 C®, etc. Halfway between the two side-walls the number 
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of calories passing through unit area changes with time according to the ex- 
pression 


x = n 
—4 — xt —16 — xt — 36 — xt 
cy=20(1-¢ ye Fe Pca 


Here the flow has already reached the fraction 0.63 of its final value at a time 
when practically no heat has started to flow at the cold end. Each section 
has its own law governing the velocity with which the steady state is attained, 
because the parts near the hot plate have to transfer energy to the colder parts, 
whereas, toward the cold side, where this task is easier, it takes an appreciable 
time for the heat to arrive. From the moment, therefore, at which in the ideal 
case one side of the slab is suddenly brought to the desired temperature, until 
the time when the temperature gradient has become constant throughout the 
thickness, the shape of the temperature waves penetrating the plate varies 
from point to point so that no definite velocity of propagation can be assigned 
to the wave fronts. Moreover, the equation is not meant to apply to minute 
intervals of time or space, such as the transfer of heat from molecule to mole- 
cule, because this exchange is not governed by k, c and s only. 

The last column of Table I gives the results for the velocity of heat trans- 
mission between an upper hot and a lower cold layer of air (k, 0.000055; 
c, 0.237; s, 0.00129; x, 0.179). On account of the low density the steady state 
is rapidly established despite the low heat conductivity k. The figures would 
not apply if a layer placed between the hot and cold “plates” were considered 
on account of the finite contact resistance to heat flow between the solid walls 
and the neighboring air layers. Apart from radiation and convection the 
quantity of heat flowing in the steady state from 1 sq. cm. of a horizontal 
plate, kept at temperature 1, into the air layers below it where the temperature 
iS Ve, is given in cal. per sq. cm. per sec., by the formula 

C= caus(04 —2») =0.000077-V/ 0, —05 (21m). 
The fourth root changes only from 1.8 to 3 when the temperature difference 
(v,;—v2) varies from 10.5 to 81° C., the range over which the formula has been 
verified (4, p. 54). On going from the solid wall into the air layers, the tem- 
perature drops fairly rapidly at first and then assumes the value which it has 
in the layers farther away from the surface. The thickness of these boundary 
layers is of the order of 0.4 cm. for air at ordinary pressures. For the range 
of temperature differences tested the heat passing from the plate into the air 
has nearly the same value as if the temperature differences were established 
between air layers 0.3 to 0.4 cm. thick, but the gradient is not linear and the 
change with time will follow a different law. On the other hand, when a slab 
of material possessing the heat conductivity k; is separated from the hot plate 
by an air space several mm. wide, the coefficient of heat conductivity k of the 
assembly is given by the expression 
1 3 ks 
Ca ; 
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The higher the value of k; the larger will be the change in &; but this change is 
negligible for very good insulators. 

As it will take nine times longer, when one plate is three times as thick as 
the other, to get the same fraction C,/C° of the steady current C®, 15 to 30 min. 
should be ample for thicknesses below 1 cm. to arrive at the steady state. 
If C, is measured at a few intervals with V maintained constant automatically, 
the value of k can be checked and the diffusivity determined at the same time. 


: Slight Annealing or Quenching of a Slab 


The same functions appear when the case is treated in which a thin slab at 
uniform temperature is suddenly plunged into a bath maintained at V° C. 
The first term in the brackets of the general formula is then the one that 
remains 


n=l 


a 
: * tee, 
oi Dd € "e* cnn Fe [509 sin n>- 2 ds, 
o 

or 

av/-2s it 

v= —| € sine —+}6 ow © de Qam + ccs ’ 

b b b 

and 


dv k on 3 =o xt z 
C=—-k—=+4V—[€ YY cosxe—+E * cop3x—-+...}, 
dz b b b 


a formula which has the same value at z and (b—z). At any distance z and 
time ¢ the flow of heat through unit area from the two boundary planes into 


‘ 4 ‘ V 
the slab is the resultant of a number of sinusoidal heat waves 4k Fj 00s 7 : 


V ‘ ‘ ; : 
Ak > 00s 3n Fete. The higher their order the more rapidly do these waves 


diminish, so that after a certain interval, ¢, only the longest wave is responsible 
for the heat transfer. The series occurring in the formula is once more a 
theta function: 


V - Bw 
Catk— Os s/b,q }, with€ @ =, 


One Boundary Plane at a Fluctuating Temperature 


When, in the “‘plate’’ method of measuring heat conductivity, one side of 
the slab is brought into close contact with a surface heated by alternating 
current, a steady state will be reached, in the course of time, such that the heat 
flowing into the plate through the hot wall during one half-cycle is equal to 
the heat flowing back during the next half, and it is of interest to ascertain 
how the heat flow varies with the thickness of the plate. On the assumption 
that the temperature at the surface z=b changes according to the law 
Vcos 2nft or Vcos pt, f being the number of cycles per second, the tempera- 
ture distribution is given by the relation (1, p. 212) in which up = Vp/2x 
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vit eee io 4 tin w (lis dl 
21 sinu(1+4)d sin wp (1—4)d 
27 


2V 2 n - 2 ut 
i = (—1)"sin n—~ ¢ ss ‘ ’ 
on y 8R+S 
The flow of heat—kdv/dz through unit area of the section at z and time ¢ 


is therefore represented by the formula 


__ whVf (144) cos w(1+i)s i, (1—i) cos m(1—ds yy] _ 
P | =u wees 


xt 

2EY S(—1)"c00 ws 8S aor c= 
e + as ns? 

For the more common insulating materials x varies from 5X107* to 5X107%, 

and f may vary from 107° cycles per sec. that is, about one cycle per day, to 

100 or more cycles per sec. As 0} will usually amount to a few cm. at least 

the fraction in the second term is nearly unity for very slow changes; for 


frequencies above one cycle per sec. it : no more than a negligible fraction. 
The whole sum is thus smaller than 2k — ( 0(2/b,q) —1 ). and will vanish 


with time at about the same rate as rt steady state is established in the 
case of the “plate’’ method. After about 10 min. only the first term con- 
tributes to the transfer of heat in the slab. In this stage the flow has become 
simply periodic so that the exchange of heat, H, per half-cycle is readily com- 
puted for the two boundary planes by integrating over the interval 0 to z/p 
and z/p to 2 z/p. 
ee! (f2 cos w(1+i)2 iy _ (1-4) cos w (1a) cv) "9 
2p sin w (1+7)d sin uw (1—1)b a 


This gives for the amount of heat entering or leaving each half-cycle at the 
side of the variable temperature the following expression 
yb sin h2yub—sin 2yubd 
m= £3 ? cos h2ub—cos 2 ub : 
As the thickness of the slab considered becomes greater, the exchange of 
heat tends towards the constant value + k V V1/zxf in agreement with the 
formula given by the simpler theory developed for this special case (Table IT). 


TABLE Il 
CHANGE OF HEAT TRANSFER PER HALF-CYCLE WITH b 


| 


aS aescccasc-asoc-seessa-aee-ee 
«/2 | 32/4 | 7 | Sait | 3%/2 rit | 2a 


0.690 0.75 0.92 | 1.0 1.02 1.01 1.00 


oor 
F ay lat cima 
1 
—0.128 | —0.134 M |-0. 34 |-0 398 \-°. 36 |-0.13 —0.08 
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On the other hand the amount of heat flowing during each half-cycle into, or 
out from, the cold plane is given by the equation 

cos ub sin hub—sin ub cos hub . 

(cos ub sin hub)?+ (sin ub cos hub)? 

It has its highest value for small thicknesses. The difference, finally, 
between H, and Hp represents the capacity of storing heat during one half- 
cycle and releasing it during the next half. 

Other periodic changes may be treated by combining a series of waves. 
The solution as given applies only to those cases in which the temperature 
of the boundary planes themselves undergoes fluctuations; when instead, the 
temperature of the surrounding air varies up and down, the contact or surface 
resistance which is offered to the heat flow by the surface of discontinuity 
must be taken into account. Its effect on heat transference through the slab 
will be most marked in the case of good conductors. The expression giving 
the heat flow also represents the distortion which an alternating electric current 
suffers in passing through a long transmission line or cable grounded at the 
farther end. It is necessary only to replace x by 1/ric:, where 7; is the 
resistance, c; the capacity per unit length. 


H, =+2—kV 
p 
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A CONVENIENT MECHANICAL MEANS OF WINDING 
QUARTZ SPIRALS' 


By J. S. Tapp? 


Abstract 


A mechanical device has been constructed for the automatic winding of quartz 
spirals with any desired pitch. The resultant spirals are compact, accurate and 
uniform, and are produced with a minimum of personal attention. A study was 
made ad the relationship between the sensitivity and the diameter of the fibre 
employed. 

he machine itself is inexpensive and involves no complicated parts. About 
50 successful spirals have been prepared in the manner recorded. 


Preparation of Fibres 


The most essential prerequisite was a uniform quartz fibre of the proper 
thickness. The quartz could best be drawn out by a falling weight acting 
through a system of strings and pulleys to give a horizontal pull. A foot- 
operated trigger held the weight in check and when released simultaneously 
removed the flame from the quartz. Fibres about 9 or 10 ft. long were suc- 
cessfully produced in this way and a little practice soon enabled the operator to 
predetermine, with considerable accuracy, the resultant diameter of the fibre 
by an adjustment of the amount of quartz heated. Since only 4 or S ft. of 
quartz fibre was necessary for a satisfactory spiral the unsuitable portions of 
the original 9 or 10 ft. could be discarded and in this way a sufficient length of 
the proper size could almost invariably be obtained. 


The Winding 


This operation, when performed by hand, required about 45 or 50 min. of 
very close attention and, even at best, the spiral obtained was unevenly spaced. 

To dispense with this tedious process a machine was designed and constructed 
from ‘‘Erector’’ parts to perform the same task more easily and to produce 
uniform coils. The device automatically revolved the rod around which the 
fibre was being wound and at the same time advanced the rod horizontally so 
as to give the resultant spiral a uniform and adjustable pitch. Furthermore, 
there was maintained on the fibre a friction tension which could be conveniently 
altered during the course of the winding if desired. 

The essential parts of the machine are shown in the accompanying sketch. 

The power was supplied by a small six-volt d.c. high-speed electric motor, A, 
which revolved at about 3200 r.p.m. A pinion, B, on the motor shaft engaged 
the flat gear, C, which in turn through a pulley and belt caused the shaft, E, 
to revolve with a speed reduction of 9 to 1 relative to the motor. This shaft, 
by means of a worm gear, H, further reduced the speed (25 times) and conveyed 
the power to shaft L which, through a pinion and large flat gear, caused the 
chuck, R, to revolve at about 4r.p.m. A 7 by } in. carbon rod, S, was held 


1 Manuscript received April 4, 1932. 
C Contributton from the Physical Chemistry Laboratory, McGill University, Montreal, 
‘anada. 


* Demonstrator in Physical Chemistry, McGill University. 
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in the chuck and kept in line by a free bearing, 7, at the oppositeend. Directly 
behind the large gear and turning with the chuck shaft was a three-inch flat steel 
disk, F, one face of which was covered with a sheet of thin cardboard. Against 
this disk was pressed a three-inch wheel with a narrow friction rim, D, mounted 
on a shaft at right angles to the one previously mentioned. By means of a 
pinion and crown gear (3-1) the movement of this shaft was communicated 





Fic. 1. Diagram of apparatus built for the winding of the quarts spirals. 
. (Length, about 18 in. ; height, about 10 1n.) 


to a vertical shaft, Z, extending down through the platform of the machine. 
This shaft had a worm gear attached at its lower end which meshed with a flat 
gear, Y, (25-1), mounted on a short shaft placed at right angles to the longi- 
tudinal axis of the machine. The lower part of this same gear meshed with a 
chain which was securely fastened level with the base and parallel with the 
longitudinal axis of the machine. As the chuck revolved, the system of gears 
just described caused the gear in contact with the chain to advance along its 
length at a speed determined by the radial position of the friction drive, and 
thereby carry along the whole machine at a uniform but adjustable rate of 
speed. (U; and U; are steel skids which support the machine and slide in a 
grooved track W.) 

In Fig. 2 is shown the apparatus used for heating the fibre and applying 
tension. A isan air-blast gas burner, B the revolving carbon rod on which the 
spiral is wound, C, and C, are two plates with narrow vertical slits through 
which the fibre H passes. D is the friction plate bearing against a similar one 
beneath, and E£ is the tension adjusting screw. WN is a pivot around which the 
upper part of the tension plate may be revolved to allow of convenient threading 
-of the fibre. R is another pivot which permits the elevation or partial rotation 
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of the entire guide so as to assure correct alignment at all times. /Y is a firm 
upright which is securely fastened to the wooden base and supports the guide 
and tension mechanism. 

The heat required to soften the quartz thread was supplied by the air-blast 
gas burner mounted on a retort stand and held in position by a clamp. The 
exact position and intensity of this flame required a great deal of study and trial 
before it was correctly placed. In the first place it must be of sufficient intensity 
to soften, and yet not appreciably weaken, the fibre, for in the latter case the 

tension would cause the 

fibre to thin out and snap. 

A very hot flame directed 

at not quite the right 

angle was also useless, be- 

cause once the quartz fibre 

touched the carbon rod 

no amount of heat would 

then cause it to become 

permanently bent. The 

best results were obtained 

when the peak of the blue 

Fic. 2. Diagram of the tension and guide device, inner cone of the flame 
mounted on the same base with the winding mechanism just touched the fibre about 


and placed at right angles to the plane of Fig. 1. ina . 
position opposite the carbon rod (S). 2 mm. in advance of the 
carbon rod, and when the 


axis of the flame, if produced, would pass down the centre of the incoming 
fibre. Many other arrangements were tried but none met with as good 
success as the one just described. 

Since the amount of extension caused by any given force is directly propor- 
tional to the number of turns in the spiral and also to the diameter of the spiral, 
it would be advantageous to have as many turns as possible, and for the sake 
of compactness, as close together as possible. Increasing the diameter of the 
spiral proportionately reduces the number of turns that can be obtained from a 
specified length of quartz. In the author's opinion, the spiral with the small 
diameter and more numerous turns has an advantage in all cases, and certainly 
has in investigations carried out under high pressures where smallness of the 
container is a safety factor. The machine was capable of adjustment to give 
spirals with from 20 to 50 turns per inch as desired. The optimum was found 
to be about 42 or 43, any more than that resulted in occasional overlaps 
which were ruinous. 

The spiral when finished was slipped off the carbon rod and neat rings formed 
at each end. 


Calibration 


The process of calibration consisted merely of suspending the spiral from 
some solid support and measuring with a cathetometer the length between the 
upper tip of the top ring and the lower tip of the bottom ring. This length 
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for any spiral was termed its ‘‘normal length’’. Small calibrated metal weights 
were then suspended from the bottom and the length measured as before. 


Sensitivity 


The weight in grams divided by the extension in mm. was termed the 
“sensitivity”. The sensitivity, when multiplied by the fraction of a millimeter 
to which the cathetometer was capable of measuring accurately (; mm. in 
this case) gave the “‘limit of detection” for the spiral in question. 

The “maximum load’”’ was that weight which caused an extension of two and 
one-half times the normal length. 

Several experiments were carried out with spirals made of fibres of different 
diameters and it was found that, for a spiral of 80-90 turns, the relationship 
shown in Table I held approximately. All spirals were of the same diameter, 
about # in. 

TABLE I 


RELATIONSHIP BETWEEN DIAMETER OF FIBRE AND SENSITIVITY OF SPIRAL 
(DIAMETER OF SPIRAL, ABOUT } IN.) 


Diameter of fibre, mm. 0.13—0.15 | 0.10 —0.11 0.09 —0.10 0.08 —0.09 
Sensitivity, gm. per mm. <0.03 | 0.011—0.018 | 0.0072—0.0080 | 0.0040—0.0060 
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THE ELECTRICAL CONDUCTIVITY OF AQUEOUS 
SOLUTIONS OF HYDROGEN SULPHIDE AND 
THE STATE OF THE DISSOLVED GAS' 


By R. H. Wricut? anp O. MAass? 


Abstract 


This paper is the last of three dealing with the equilibria between hydrogen 
sulphide and water, and is a continuation of a series of researches designed to 
investigate the equilibria existing in gaseous aqueous systems. 

Electrical conductivities of aqueous solutions of hydrogen sulphide have been 
measured between 5° and 60°C. The results lead to the conclusion that hydrogen 
sulphide resembles certain other gaseous solutes and forms with water a complex 
which undergoes electrolytic dissociation. The constant & of the Ostwald dilution 
law therefore appears to be an apparent, rather than a real, dissociation constant. 

The complete analysis of the equilibria involved and the evaluation of the 
constants must depend on accurate measurements at low concentrations, which 
are not yet completed. 


Introduction 


The state of hydrogen sulphide dissolved in water has hitherto been known 
only from the scattered observations of a number of observers. It has been 
the object of the work described in this and two preceding papers (10, 11) to 
secure an homogeneous body of information from which the nature of aqueous 
solutions of hydrogen sulphide could be deduced. To this end, the partition of 
hydrogen sulphide between the vapor and aqueous phases has been deter- 
mined and the conductivities of the solutions have been measured. This paper 
is devoted to the latter of these measurements and the preliminary theoretical 
consideration of the problem as a whole. 

The object of the partition experiments was to determine the solubility curves 
in such a way as to eliminate all sources of deviation from the ideal laws not 
arising in the liquid itself. For instance, an apparent departure from Henry’s 
law (regarded as a particular case of the partition law) might easily be caused 
by failure of the hydrogen sulphide in the vapor to obey the simple gas law*. 
Any departure, therefore, from the partition law after allowing for these 
extraneous effects may tentatively be ascribed to conditions in the liquid such 
as association, dissociation, hydration, etc. 

The electrolytic dissociation of hydrogen sulphide is known to take place in 
two steps but, as it has been well shown (2, 5, 6) that the constant of secondary 
dissociation is vanishingly small, in this work hydrogen sulphide will be treated 
as a monobasic acid. 

*That allowance for this factor is indeed necessary 1s shown by the results in the previous fo- 
(11), in Fig. 5 of which the solid lines show the actual variation of the concentration of the solutton 


with concentration of hydrogen sulphide in the vapor, whereas the dotted lines represent the curves 
obtained when the vapor was treated as an ideal gas. 


‘ Manuscript received A pril 8, 1932. 
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The experimental data for the electrical conductivity submitted in this paper 
are the first that have been published having any claim to comprehensiveness 
or precision. 

Determination of Electrical Conductivity 


Experimental Method 

Previous attempts to measure the electrical conductivity of hydrogen sul- 
phide solutions have been hampered by experimental difficulties arising from 
polarization of unplatin- 
ized electrodes and occlu- 
sion of hydrogen sulphide 
by platinum black. Control 
and measurement of the 
concentration have also 
been difficult (1, 9). 

By adapting the sealed 
cell procedure described in 
the preceding paper (11), 
the latter of these difficul- 
ties has been overcome, 
and the problem of polar- 
ization has been avoided 7 M 


altogether by using the ts tee ia = ' 
: IG. 1. 4agram showtng arrangement of apparatus 

static method of conduc- for the electrical conductivity determinations. 

tivity measurement. This 


method is not new but it has lately been improved in these laboratories and 
elsewhere and seems quite satisfactory. 

The principle of the method is very simple. If two resistances R; and R, 
be connected in series and a current passed through the circuit, then if E, is 
the potential across R,, and E; the potential across Re, 

Ex _R 
ER 

Fig. 1 shows the application of this principle to the measurement of electrical 
conductivities. The conductivity cell, C, was connected in series with a 
resistance, R, milliammeter, M, battery of dry cells, B, and a switch, S. The 
cell was provided with two exploring or ‘‘secondary"’ electrodes which lay 
between the primary ones. 

By means of the double pole, double throw switch, D, the secondary elec- 
trodes or, alternatively, the resistance, R, could be connected to the quadrants 
of a Dolezalek electrometer, E. Since the deflection of the electrometer was 
related to the potential difference of opposite pairs of quadrants, then, from the 
ratio of the deflections the ratio of potential fall between the exploring electrodes 
to the potential fall across the known resistance could be obtained. Since the 
current was the sarne through both, this ratio was also the ratio of the re- 
sistances. The absolute values of the potentials and of the primary current 
were within limits immaterial, so that the resistivity of the solution was found 
in terms of two electrometer deflections and a known resistance. 


E 
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Since, however, the electrometer deflections were for various reasons not 
strictly proportional to the applied potentials, the electrometer was first 
calibrated by connecting a number of combinations of standard cells to the 
quadrants and plotting a curve of potential against deflection. 

The usual precautions of insulation, shielding, etc., were taken in setting 
up the apparatus. * 

The conductivity cell used is shown in Fig. 2. It was made of Pyrex glass 
and the electrodes were of thin platinum ribbon, unplatinized and adhering 
closely to solid glass sup- 
ports as shown. The sec- 
ondary electrodes, B and 
C, had a very much smaller 
surface than the primary 
and consisted simply of 
small platinum tips pro- 
jecting into the side arms, 
the entire side arm con- 
stituting the exploring 
electrode. 

The cell was mounted 
on a mechanical rocker 
placed in a_ thermostat. 
The cell contents could be 
thoroughly stirred by rock- 
ing the cell back and forth 

Fic. 2. The conductivity cell. between readings. The elec- 

trometer connections to the 

mercury wells were made through flexible wires and were not removed during 
the rocking. 

A cell constant was determined using 0.02 N and 0.01 N KCI, the specific 
conductivities of which were taken as 0.002765 and 0.001413 respectively, at 
25° C. (7, p. 213). The water correction was obtained from a blank run. 

The procedure in filling the cell was similar to that described in the previous 
paper (11). The water was not, however, run directly into the conductivity 
cell from the weight pipette, but into a bulb temporarily sealed to the apparatus. 
In this bulb it was freed from dissolved air by freezing and melting im vacuo and 
was then distilled into EZ, Fig. 2. The hydrogen sulphide was purified and 
measured as previously described, and condensed in F with liquid air while 
the water in E was frozen with solid CO). The cell was then sealed off at G 
and the conductivity of the resulting solution measured at various temperatures. 
At the conclusion of the measurements the cell was opened and its volume 
found by filling with distilled water and weighing. 





*It ts usual to ground one pair of quadrants, but it was found in this work that better results 
were oblained by grounding a point in the primary circuit instead. 
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Results 


The concentration of the solution at each temperature was found as shown 
in the following specimen calculation. As the partition coefficient of hydrogen 
sulphide between liquid and vapor varies with the pressure, an approximate 
calculation was first made showing roughly the pressure in the cell. The 
appropriate value of the distribution coefficient was am selected from the 
data of the previous paper (11, Fig. 6). 


Specimen Calculation 

Temperature, 20° C.; wt. of H.S, 0.7439 gm.; wt. of H,O, 42.38 gm.; vol. of 
water, 42.49 cc. ; approximate pressure, 2100 mm. ; approx. wt. of H,S in solution, 
0.460 gm.; increment of volume due to dissolved gas (from mixture rule), 
0.54 cc.; volume of solution, 43.03 cc.; volume of the cell, 113.08 cc.; volume 
of the vapor, 70.05 cc.; partition coefficient D (20° C. and 2100 mm.), 2.67. 

D X (vol. of soln.) X (total wt. H:S) 

(vol. of vapor) + D X (vol. of soln.) 

2.67 X 43.03 XK 0.7439 

70.05 + 2.67 X 43.03 
= 0.462 gm. 


0.462 x 1000 
34.08 X 43.03 


= 0.0316. 


Weight of HS in solution 


Molarity of solution, 


From the scale readings, electrometer calibration curve and cell constant, 
the specific conductivity was obtained at a number of temperatures, and values 
at other temperatures were interpolated. To obtain the degree of dissociation, 
the limiting equivalent conductivity at infinite dilution was calculated from 
the data of Jellinek and Czerwinski (4) who measured the conductivity of 
NaSH solutions at 0°, 18° and 25°C. The ion conductances of Nat and H* 
were taken from the International Critical Tables (3). The limiting con- 
ductance of hydrogen sulphide at other temperatures was obtained by inter- 
polation. These conductances are given in Table I. 


TABLE I 

LIMITING EQUIVALENT CONDUCTANCE OF HYDROGEN SULPHIDE 

Temp.,°C.| 0 5 10 15 18 20 
A 269 | 299 | 328 | 356 | 372 | 385 


25 30 40 50 60 
414 | 442 | 499 | 556 | 613 





The constant, k, of the Ostwald dilution law (i.e., the apparent dissociation 
constant) was calculated from the formula, 


AeX* 
u(t) 
ba __ \4e/_ 
Ay 
Ao 


1- 
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where, M is the molarity of the solution, A, the equivalent conductivity at 
concentration M, and A, the equivalent conductivity at infinite dilution. 

The results are summarized for each temperature in Table II, and in Table III 
values of the apparent dissociation constant are collected and compared with 
those of other authors. 









































TABLE II 
EXPERIMENTAL RESULTS 
| | 
Total Partial . | Specific Ay 
pressure | pressure Molarity | cond. X 10° ‘ x 10° kx10* 
"cin as Yi, ‘ 2 SR 
O7 601 (0.141 2.4 ~«|~«0.578 ) 4.72 
903 | 897 . 208 | aoe 479 | 4.77 
1181 _ 1175 me ao? oe Sy 415 | 4.65 
At 10°C. 
684 eae Ee (re )3=— rr 
1011 | 1002 ae ok. | 547 | 5.98 
1321 | 1312 | .259 | 4.12 | 485 | 6.09 
me | feet ee he ok 406 | 5.64 
At 15°C. 
te ae | a ae | ee 
1117 1104 .190 | 4.28 631 7.59 
1459 1446 . 248 4.92 557 | 7.70 
1940 1927 | 329 5.45 | 466 72 
At 20°C. 
a sis j|t«iéi 4.04 | 841 | 8.75 
1224 1207 =| 182 4.98 | 710 2=— | 9.19 
1611 1594 240 Pe i 619 | 9.20 
sigs enindinnsi eee ich S lO kncdcit RMaliaille Rs 
At 25°C. 
907 883 118 4.64 | .948 | 10.6 
1332 1308 174 | 5.72 | 794 11.0 
1741 1717 | 228 | 68s | 694 | 11.0 
2308 | 22h | 302733 | 88610 
At 30°C. 
982 950 tn aan 1.05 | 12.4 
436020] 140g 167 | (6.48 s79 «=| s(12.9 
1879 1847 218 te 767 | 12.9 
2480 | 2448 =| «289 | Oe | 650 12.2 
At 40°C. 
Sn he | | nd 1.24 | 16.2 
1639 1584 153 | = i ee ee 
2147 2092 202 | 9.85 | 910 | 16.7 
2829 2774 267 10.3 — a). 
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TABLE Ill 
SUMMARY OF DILUTION LAW CONSTANTS 





| Ff | 








| | | | | 
Temp., °C. | S | 10 | 15 | 20 | 25 | 30 | 40 | 50 | 60 
k, mean value X 108 | 4.71 | 5.74 | 7.47 | 8.96 | 10.8 | 12.6 | 16.4 | 20.3 | 23.9 
Values from the literature 
Nl 
, Observer | Jellinek and Auerbach (1) Walker and 
Czerwinski (4) Cormack (9) 
Temp.,°C. | 0 | 18 18 
k 1 x 10-8 9.1 x 10-8 | 5.7 x 10-8 
Discussion 


The most striking feature of the experimental results is the regular and 
marked increase in the dilution law constant, k, with rise in temperature. 
_Applying the equation 

dink —@q 


dT RT?’ 
+q will signify the heat evolved in the reaction 
HS —»> H* +SH +4. 
The ordinary method of integration can be applied using the formula, 
k 1 1 
we 8- (a > H)- 

Taking values of k from Table III, it is found that integrating from 5° to 
25°C. gives +q=—6800 cal., whereas integrating from 30° to 60°C. gives 
+q= —4300 cal. 

So great a change in g can hardly be explained by Kirchoff’s theorem in 
terms of changes in the heat capacities of the ions or the dissolved hydrogen 
sulphide, and another explanation must be sought. 

By assuming that when hydrogen sulphide dissolves in water, the solvent 
and solute unite to form a complex which then undergoes electrolytic dissocia- 
tion, the above result becomes intelligible. The quantity g will then depend 
on two independent quantities, Q, the heat evolved in forming the complex 
from hydrogen sulphide and water, and Q2 the heat evolved in the ionization 
of this complex. If Q: has a negative value, i.e., if the dissociation of the 
complex is attended by the absorption of heat (which is to be expected from 
the similarity of hydrogen sulphide and water, an increase in the extent of 
ionization would be expected with rise in temperature. This tendency would, 
however, be opposed by the decrease in the dielectric constant of water, and in 
any case the degree of ionization seldom varies extensively with temperature. 

The observed increase in k may also be explained in another way. Supposing 
Q, to have a negative value, i.e., that the H»S-H,O complex is endothermic, 
then a rise in temperature would favor its formation and would therefore produce 
an increase in the number of ions over and above any increase due to a change 
in the ionization constant of the complex. 
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Furthermore, it is evident that —g is the net heat evolved when (a) the ions 
combine to form the complex and (b) when a part of the complex then dis- 
sociates into free hydrogen sulphide and water. Assuming an endothermic 
complex, the process (b) will be relatively more important at low than at high 
temperatures in so far as it contributes to the net heat evolution g. The 
marked decrease in g at higher temperatures is therefore readily explained by 
the hypothesis of an endothermic complex. 

The system may now be provisionally formulated in the following way 
(secondary dissociation being neglected). 


H,S (in vapor) “= HS (in solution). 
HS (in soln.) + »H,O “= H,S.nH:0O + Q:. 
H:S.nH.O “> H*t + SH~.nH:O + Qs. 


Corresponding to each equilibrium there will be an equation: 


[HS] soln. [H,S] [H,0]” [H,S.nH,0] 
[Hs] vap.  ” (H.S0H,0) ~~” {H*) [sH-.»Ho] ~ 


It is evident that the system is now analogous to the systems treated by 
Maass and Morgan (8), viz., sulphur dioxide and water, carbon dioxide and 
water, and ammonia and water. There is however an important difference. 
In all the systems treated by these authors, the formula of the complex 
(H,SO;, H2CO;, and NH,OH) could be inferred from the composition of the 
salts formed, whereas there is no such clue to the value to be assigned to in 
the equations above. 

Maass and Morgan found that measurements at high concentrations were 
necessary for the evaluation of K; and K;: in the systems examined by them. 
From a consideration of the present problem, it appears that measurements 
with very low concentrations are also required when 7 has also to be determined. 
Since these have not yet been made, this paper can include only a qualitative 
consideration of the system. 

The assumption of a complex such as has been described will also account, 
in a qualitative way, for the form of the vapor pressure curves of the preceding 
paper (11, Fig. 5). From the equation, 


(H,S] [H.0]” 

ee 
it is evident that as [HS] increases (at constant temperature) [H:S.nH;O] will 
also increase, but relatively more slowly owing to the diminution in [H,O]. 
The net result of [H:S] increasing more rapidly than [H:S.nH,O] will be a 
curvature in the solubility curves of the kind actually displayed by the iso- 
therms of the previous paper (11). 


Ki 


References 


1. AvgrBacn, F. Z. physik. Chem. 49: 217-223. 1904. 
2. Aumfras, M. Compt. rend. 186: 1724-1726. 1928. 
3. INTERNATIONAL CRITICAL TABLES, v. 6, MCGraw-Hiti. 1929. 





1 
1 


KB OOaOnN AM 


CONDUCTIVITY OF AQUEOUS HYDROGEN SULPHIDE 


. JELLINEK, K. AND CZERWINSKI, J. Z. physik. Chem. 102: 438-479. 1922. 

. Knox, J. Z. Elektrochem. 12: 477-481. 1906. 

. Knox, J. Trans. Faraday Soc. 4: 29-50. 1908. 

. LANDOLT-BORNSTEIN. Physikalisch-chemische Tabellen. v. 2. Springer, Berlin. 1923. 
. Morcan, O. M. and Maass, O. Can. J. Research, 5: 162-199. 1931. 

. WALKER, J. and Cormack, W. J. Chem. Soc. 77: 5-21. 1900. 

. Wricut, R. H. and Maass, O. Can. J. Research, 5: 442-447. 1931. 

. Wricat, R. H. and Maass, O. Can. J. Research, 6: 94-101. 1932. 





AN EQUATION OF STATE FOR GASES AT 
LOW DENSITIES' 


By D. LEB. Cooper? anp O. Maass?* 


Abstract 


An equation of state for gases at low densities is developed, using a new function 
for the change in viscosity with temperature, also developed herein. 
The gas law equation takes the form 
RTW1+KT) _ 


PV?+a—RTV— i+KT 0 


or V(i+KT)(PV—RT)=AT—a where a and 3b are constants corresponding to 


those of the Van der Waals’ equation, and K is a constant derived from the 
proposed viscosity function which is, for carbon dioxide, 


g=v T(1+KT) 
where XK is a constant and 4 is the viscosity at ap absolute temperature T. 

In the case of carbon dioxide the equation was found to follow density results 
with an accuracy of from 0.01% to within experimental limits, and the viscosity 
oo was found to agree with Sutherland's (10) results between —78.5 and 
2 5 

Comparisons with several other equations of state are made. These show 
that the new equation is probably more accurate than any other. 

An expanded form of the new equation, namely: 


aid RTby ia RTbo \2 
M ~14( I+KT )P +2( +K7 ) P + etc. 
RT? RT? 
permits calculations of the slopes of isothermals for any temperature. Com- 
parisons are made with experimental data. 
The expanded form of the equation may be solved for K, giving the expression: 
(i—%) (A+K) — AA 
(0-83) (42+K) AAs 


2. wie & es 7 - 
where 8 T and A=a o+K and ¢ = Rb, and since the equation enables the 


calculation of the molecular radius 7, the viscosity may be calculated for any 
temperature and pressure over which the equation holds. 


Introduction 


The authors (5) have published the results of measurements of the density 
of carbon dioxide. The accuracy of the results was of the order of 0.01%, 
approximately 10 times that of previously existing data. The object of this 
communication is to apply these previously published results to the examination 
of some equations of state including a new one proposed herein. 

The application of the new equation is restricted to gases at low densities. 
The restriction is intentional and succeeded a consideration of the possibilities 
of the usefulness of equations of state for deductions concerning molecular 
phenomena. These considerations are outlined below. 

An infinitely dilute gas is usually considered ‘‘ideal’’. More correctly, it 
approaches the ideal as its density decreases. Increasingly large aberrations 
follow increasing density. Representation of these aberrations by correction 

+ Manuscript received A pril 19, 1932. ‘ on 
me from the Department of Physical Chemistry, McGill University, Montreal, 


2 Holder, at the time, of a fellowship under the National Research Council of Canada. 
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terms applied to the ideal gas law equation presupposes at least partly analogous 
conditions between the actual and the ideal. Representation of high pressure 
data is complicated by the fact that neither is the term ‘‘molecular volume”’ 
unique, nor ‘mean free path”’ significant. As a consequence it appears to be 
a reasonable supposition that more pregnant deductions can be made by the 
use of low density data and equations holding in that region in which partially 
analogous conditions exist. 

The proposed equation is developed by the use of a Van der Waals’ (12) 
molecule. Possibly that developed by Lennard-Jones (8) has the greatest 
significance compatible with tractability, but with a restriction to the low 
density region and assumption that the effective molecular radius does not 
vary within the second order of small quantities, and furthermore that at low 
densities the attractive field is not appreciably affected by pressure, the use of a 
simpler Van der Waals’ molecule is not only justified but preferable. This is 
more particularly true in the treatment of gases with large molecular fields 
and relatively soft molecules. 

The degree of conformity of the equation based on these postulates deter- 
mines whether or not the above assumptions are justified. 


Mathematical Development of Equation 


Maass and Mennie (9) have shown that the pressure of an ideal gas may be 
represented by the equation 


(P+ 7) V= (1+ 7) RT, (1.0) 
where r is the molecular radius, / the mean free path, A a Van der Waals’ 
constant of attraction, and P, V, R and 7 have their usual significances. 
If, as postulated above, r and A remain effectively constant, the variation in 
the deviation (P V-RT) results from changes in /. 

The variation in / may be represented as a function of the temperature. 
Sutherland’s (11) expression incorporated in equation (1.0) yielded the 
following : 
PV!—RTV+A -RTB(1+ 5) =0, (1.1) 
x 


where 


_ 8V2 «FN 
1+ 


. x a (1.2) 
273 


PV?—RTV+A —RTB'=0. (1.3) 


Maass and Mennie (9) point out that this equation may be transformed into 
one identical with that of Van der Waals’ by a first order approximation, and 
that in agreement with the fact that Sutherland's function fails at low tem- 
peratures, Equation (1.3) fails at temperatures below 0° C. 

It was believed that an equation of this type was capable of further develop- 
ment provided a rigid function for variations in / with temperature could be 
found. The following is offered as a solution. 
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The mean free path / is related to the viscosity by an equation such that 
n=Kymnzl, (2.0) 


where K, is an unknown constant. Or, more accurately, since 7 changes more 


rapidly than /T 
” = Kymnxif( 7. (2.1) 


But yeVT, 
therefore n=Kymnv/T f(T), (2.2) 
i ala aes 
4/2 rn’ 

where y is some constant, and as /=4 when 7=% (assuming as before 
that the total change in 7 is caused by a change in /) therefore 

_ Kimn/ THT) 

4V 2 arn 

whence n=KsVT f(T), (2.4) 


K 
where K;= os =const., (2.4.1) 


and the viscosity at0°C.is 4h 


’ (2.3) 


=K;f(T). (2.5) 


” 
or JT 


The data of Sutherland and Maass (10) was used to plot VF against T 


for carbon dioxide and f (7) was found to be a straight line, whence 
ie 
VT 1+XT, (2.6) 


; 1+KT 
and we may write l=h 4 eT ’ (2.7) 


whence substitution in (1.3) gives 


i: 2r(1-+K Te) 
(?+-7) V er | 1+ ten |: (3.0) 


where a is a Van der Waals’ constant of attraction differing numerically from 
A (Equation 1.3). 
¢ 
But h 4/2 xrin’ 
Noa 8V2 xrn(1+KT) 
ee (P+ 4a) ¥ ar| 14 VI+RT) | 
RTb(1+KT») 
1+KT 
_ 82 arn 
— 


or PV?+a—RTV— =0 


and b 


Equation (3.2) may be written 
V(1+KT)(PV—RT) =hT—a+eKT, 


where bp = b(1+K To) =const, 


or Vii+KT)(PV—RT) =iT—a, 
where A=bo+aK =const. 
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Equation (3.4) is the one used for comparisons. 
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The constants K and A 


may be evaluated easily and the equation tested for rigidity without difficulty, 


as inspection shows that the 
left hand term should vary 
linearly with JT. That this 
is so is seen from Fig. 1, 
which shows that the equa- 
tion holds over the measured 
range with the desired ac- 
curacy. The values of A and 
a may be determined directly 
from the line of Fig. 1. 

The function may be eval- 
uated as follows. We define 
a volume, V, such that V 
represents the actual volume 


560 


340 


240 


$6 58 60 G2 


64 66 68 70 i 


VI+KTXPV-RT) 


Fic. 1. The calculated values of V(1+KT) 
(PV—RT) plotted as a function of T in °A. 


occupied by a theoretical molecular weight of a gas at a corresponding T and 
P. Vis then calculated from the published data from the equation, 


(A) 


For purposes of correlation these values are shown in Table I. 


TABLE I 
APPARENT MOLECULAR WEIGHT AND VOLUME 
OF A GRAM MOLE OVER A 
TEMPERATURE RANGE 


Apparent 
molecular 
weight (M,) 


44.422 
44.295 
44.232 
44.167 
44.138 


Temperature 
° Abs. 


V 
(observed) 
litres 


243.18 
273.18 
293.18 
323.18 
343.18 


19.764 
22.266 
23 .930 
26.418 
28.071 


The values of RT may be calculated 
for each temperature and those of 
1+KT from the viscosity line (Fig. 1) 
the constants of which are, for carbon 
dioxide, 1+4.09410~*7. The values 
of the left hand function lie on a 
straight line (Fig. 1) the equation for 
which is 

AT —a=1.6167 X10-*— 11.34. 
The numerical values of the results 
are shown in Table II, and the exacti- 
tude of the equation is demonstrated 


further by a direct calculation of V shown in comparison with those cal- 
culated from several other equations (Table IT). 
TABLE II 


COMPARISON OF RESULTS, V obs. AND V calc. AT VARIOUS TEMPERATURES; 
PRESSURE CONSTANT AND EQUAL TO ONE ATMOSPHERE : 


Temperature V obs., 
° Abs. litres 


243.18 
273.18 
293.18 
323.18 
343.18 


19.764 
22.266 
23 .930 
26.418 
28.071 


V;*, 
litres 


19.764 
22.266 
23 .930 
26.418 
28.071 


V2, 
litres 


Va, 
litres 


19.772 
22.265 
23 .930 
26.418 
28 .069 


19.952 
22.414 
24.055 
26.517 
28.157 


*V, =new equation; V:= Maass-Mennie equation; V;=gas law. 


74 76 
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The equation for carbon dioxide is 
V(1+4.0983 xX 10-*7) (PV —RT) =1.616 X 10-*T — 11.34, (3.5) 

or, written in a form to show minor constants, 
0.3596(1 +4.0983 x10~* T») 


PV?+11.34-—RTV—RT 
1+4.0983 x10 ~*T 


=0 (3.6) 
Discussion of the Equation 
The equation may be written in the form of an isothermal by substitution 
of equation (A) and expansion by the binomial theorem. Then, 
in RT o— RT \: 
a, 71+ ie ~14KT)P + a( +KT ) P+ ai (4.0) 
RT? _ RP 
where by = 6(1+KT7>). 
The slopes of the isothermals may be calculated from the coefficient of P. 
These are shown compared with those determined experimentally in Table III. 
A variation of 0.01% in the highest 
TABLE Ill molecular weight on any isothermal 
SLOPES OF LOW PRESSURE ISOTHERMALS results in a change of slope of 18%. 
~ The calculated values are well within 
a. . ed) | satiate ed) the limits of error. The coefficient 
—_____|______—_ of P* allows for curvature of the iso- 
thermals. It becomes effective at five 
243 .18 ° 0.413 
273 18 0.288 atmospheres to the extent of 0.01%. 
ig co Following are similar expansions of 
343.18 0.135 the equations of Maass and Mennie, 
__._ and Van der Waals, shown for com- 
Parison: 





Maass and Mennie, 
1 
Af =1+ (“ar Aan)? + “B sor) Pt. 
0 


Van der Waals’, 


A second method of treatment of the expanded form (4.0) will be discussed 
later. 


Comparison of Equations 


Table II and succeeding tables serve to compare the several equations. 
TABLE IV Table II shows the calculated volume V 
es aan i compared with those determined experiment- 
AN ISOSTERE AT V =23.930 1. ally. The values of the Maass-Mennie equa- 
tion were calculated using their method with 
constants derived from the published results 


| of the authors. The discrepancy at —30° C. 


Pobs. | P calc. 


yd 2 tones is well outside the limits of error. 


1.0000 | 0.99970 Table IV shows comparisons with Bridge- 
: oar ae man’s equation for carbon dioxide, calculated 


in the form of an isobar. At the lower tem- 
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peratures the deviations are greater than experimental errors. 


TABLE V 
COMPARISON OF THE OBSERVED PRESSURES WITH THOSE CALCULATED FROM VAN DER WAALS’ 
EQUATION USING DIFFERENT VALUES FOR THE CONSTANTS @ AND 8 














Re a en oe ae 











Temp., ° Abs. P. obs. B 
















0.166 | | 0.0712 


P. calc. 









.0000 1.0013 1.0000 1.0003 





243.18 1 1.0024 
273.18 1.0000 1.0005 1.0000 0.9999 1.0014 
293.18 (1.0000) (1.0000) (1.0000) (1.0000) 1.0009 
323.18 1.0000 0.9996 1.0000 1.0001 1.0004 
343.18 1.0000 0.9995 1.0000 1.0003 1.0004 
























Table V requires further discussion. Attempts have been made to calculate 
the constancy of Van der Waals’ a and b, written here, a and 8. Algebraically 
they are unique, and hitherto it has been impossible to detect differences 
between calculated and experimental values even with considerable differences 
in the values of the constants used. The more accurate results show that 
a 7% variation in a causes a noticeable difference at the extreme temperatures 
even when the volumes are calculated using a corresponding (. 

The values of a and # were obtained by a simultaneous solution of Van der 
Waals’ equation using the authors’ results: a=7.06; @ =0.166. Pressures cal- 
culated from these values follow the experimental results (column 4). Those 
calculated (column 6) from the best representative values of the constants 
based on the results of Amagat (1, p.109) and Andrews (2) and quoted by 
Jellinek (7, p. 632) show discrepancies varying from 0.2% at low temperatures 
to 0.03% at the highest. The figures of column 3 were calculated by substitu- 
tion of the best representative value of a and calculation with a corresponding 8. 

The results in Table V show that a change in the constants of Van der 
Waals’ equation makes it inapplicable to calculations of the highest accuracy. 

Further comparison of Equations (1.3), (3.2) and Van der Waals’ is made 
possible by the fact that the equations are similar when 

ap RboT 
RTB'=PVB— V * T+KT' 
and, with the first order approximation that PV =RT, ‘ 


RTB'=RTB— - a 


The term “© affects the results to —4— and neglecting this, 
Vv 6000 

















bo 
BinB= TTRT ° 
or at 0° C. 
Bi=g=b. 
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TABLE VI a eae 
iiaiaidadieedis ame i, Sa Rill This indicates that at all temperatures 


HIGHER PRESSURES other than 0° C. the proposed equation 
a eee ——=== demands a variable, and the others a 


Pressure, | Pressure, | Pressure constant correction factor. 


obs. (Van der W. aals'| (new equation) 
equation) | 


| Table VI shows the range over 
; which the new equation holds com- 
12.01 | 12.12 | 12,012 pared to that of Van der Waals’. The 
. = | oa a values for the latter were calculated 
20 01 | 20.74 19.945 by use of constants obtained from 
4.0 | 30 . 86 | 33.297 data of the authors. The constants 
aptioee = for the proposed equation were ob- 
tained, as explained, by a process yielding unique values, for if the value of 
K be changed 100%, i.e., be placed equal to zero, calculated and experimental 
values do not agree as is shown in Table VII. 
The forms of the equations allow TABLE VII 


another differentiation. The new equa- SHOWING EFFECT OF A CHANGE IN K ON 
AGREEMENT BETWEEN OBSERVED 


tion permits calculations of the slopes AND CALCULATED PRESSURES 

of the isothermals as shown above, "=| 
but Van der Waals’ equation has no Temp., 
like application due to the substitu- ° Abs. press. 
tion of PV by RT. 

Comparison with the equation of 
Bridgeman (4) is more difficult on 
account of the relatively large number -18 | 1.0000 1.00150 | 0.9999 
of semi-empirical constants therein. ee - ee 
Table IV indicates that while somewhat similar in form, needs equation 
fails at low temperatures and pressures. At higher temperatures and pres- 
sures the agreement is satisfactory. 

A more apparent difference between the three equations first mentioned 
may be detected by writing them in the form of isosteres, whence, for Maass 
and Mennie, and Van der Waals’ we have 

P=BT-A, 


and for the recently proposed equation, 





| Obs Calculated pressure 


K=0 ee 
| 1 


| 
Daria 

1.0000 0.99884 
| 





P=B,T+ - —a. 


TERT 


The first two demand straight isosteres, the latter allows for a curvature. 
In this respect the new equation probably represents the facts to a greater 
degree of exactness since isosteres are known to possess a curvature at higher 
pressures. No curvature was detected in the isothermals, hence no value of 
C, can be calculated. An upper limit was found to be of the order of 110~*. 


Further Applications of the Proposed Equation 
Among others the equation lends itself to two further applications, namely, 
the calculation of the molecular radius, r, and the viscosity constant K, as 
defined by Equation (2.6). A knowledge of r permits the calculation of the 
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viscosity of the gas under examination at 0° C., which value, substituted in 
Equation (2.6), permits, with a knowledge of K (also calculated from gas law 
relationships by use of the new equation), a calculation of the viscosity at any 
temperature over which the equation holds. In so far as the second application 
is concerned this is the first time that the calculation of viscosities at different 
temperatures has been possible from PVT results only. 


The Molecular Radius 
The calculation of the molecular radius may be carried out as follows: 
From equation (3.1) 


RTb(1+KT>) 
246 — , eae i 
PV?+a—RTV A+RT) 


where b=8/ 2 arn, 
therefore b=6,/2 Vo, where Vp is the total volume of the molecules; y is a 
correction factor. Using Jean’s calculation of y, we have 
_pry— V2 VAltKT) _ 
anneal) “aay 
from which r= 2.28 107° cm., the calculated viscosity at 0° C. of 0.0001353 
differing from the experimental value of 0.0001354 by 0.15%. 


The Viscosity Constant 

The second application, namely, the calculation of the viscosity constant 
K from PVT data, may be carried out as follows. 

From equation (4.0) we have 


Mo 


but ; M! = M,+uP, 
where u is a constant 


RT? 


a RT be 
ca -1+( 147) P, 


therefore 


and equating equal coefficients, 


where A= 37 R°T?, and § = Rbo. 

To solve for K, we define some quantity @ such that 0 = + 
then A=a— etx ’ (6.3) 

a iC) 
and eliminating a 4K) @,+K) AAs. (6.4) 
&(0:—0;) = 

(4+K)(%+K) 
(i—%)(&4+K) _ Ae 
(0, —63)(2+K) AAs’ 
and since 4 may be evaluated for each experimental isothermal, K may be 
calculated directly. 


Similarly ha —de, (6.4.1) 


and by division (6.5) 
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The nature of the equation is such that K may be evaluated to one significant 
figure only, its value being 5X10~*. Using this value for a calculation of b 
we have, b=0.412, from which the viscosity at 0° C. is 0.000131, compared 
with the measured value of 0.0001354, and using the value of K above, 6 at 
— 60° C. is 0.000116 against a measured value of 0.0001160 for carbon dioxide. 
Thus without the use of any results other than those used to determine the 
slopes of the isothermals, viscosities may be calculated over a long temperature 
range with an accuracy of about 2%. The nature of the viscosity temperature 
function has, however, to be assumed. 
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THE CONDENSATION OF CERTAIN y-KETONIC ESTERS 
WITH AROMATIC ALDEHYDES'! 


By C. F. H. ALLEN? AND G. F. FRAME®* 
















Abstract 


The condensation of methyl and ethyl a-phenyl-8-(para-substituted) benzoyl 
propionates with benzaldehyde and piperonal in the presence of sodium methyl- 
te, followed by acidification, has been found to produce cyclic compounds; the 
atter are shown to be lactols, six of which are described. The spontaneous ring 
closure is probably due to the highly branched chain. A mechanism for the 
reaction is proposed. 

On oxidation with chromic acid, a-diketones are formed; the latter contain the 
aryl group introduced as aldehyde. The lactols resemble triphenylcarbinol in 
several respects, forming a chloride, methyl ether, and acetate on appropriate 
treatment, all of which, on hydrolysis, regenerate the lactol. They also give 
colored solutions with concentrated sulphuric acid, from which the starting 
material is recovered on addition to water. From this work it is evident that 
in arylated y-ketonic esters the hydrogen atom alpha to the ketone carbonyl 
group is more active in alkaline aldol condensations than the hydrogen in.the 
alpha position to the carbalkoxy — group; the observation of others, that 


LO 
GY a 
the conjugated system > C==C—-C——-OR is more stable than > C=C—C—R, 
has been confirmed. 










In a previous investigation (12*) it was found that when benzaldehyde 
reacted with a-phenyl-@-benzoylpropionitrile (I) in an alkaline medium, the 
hydrogen in the position alpha to the cyano group was involved: 
C,Hs;CHOH C.H;CHO + C.H;CH 
not _d 

CsHsC—CN fe CsHsCHCH,COC,H; a> CsH;CH—C—COC,Hs 
H:COCsHs N N 

II I 
In the preparation of a-phenyl-8-benzoylpropionitrile by treatment of ben- 
zalacetophenone with an aqueous solution of potassium cyanide, a high melting 
insoluble solid is the principal product unless the hydrogen ion concentration 
is carefully adjusted; this solid has been shown to be III, formed by addition 
of the nitrile first formed to a second molecule of unsaturated ketone, and loss 


of water (8). 
C.HsCH=CHCOC,H; 

























Cs,sHs;CH———-CCOCG,Hs 
; ¢ 
+ en Cs,HsC—CH:—CC,H, + H;0 
CN N 





III 
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*In the first sentence, 3rd line, page 1353, the word ‘‘nitrile’’ ts ambiguous; it should read ; 
a-phenyl-B-benszoyl propionitrile. : 
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Thus the hydrogen atom alpha to the cyanogen radical is more active than the 
one alpha to the carbonyl group. Further examples are found in certain 
8-ketonic nitriles (1, 12) that on bromination give almost entirely bromine 
substitution products of type A. 


CsHsCHCH:COC.H,.X CsHsCHCH2;COC.H«X 
YC.H«CHCN YCsH,C—CN 


r 


A 


On account of these results it seemed desirable to investigate the ester (IV) 
corresponding to the nitrile (I) and determine which was the more active 
hydrogen. It was known that the acid (V) on bromination gave the mono- 
bromo-substitution product (VI), indicating that the hydrogen alpha to the 
ketonic carbonyl group was more reactive (14). 


C.Hs;CHCH,COC,H,X—4 ee CsHs;CHCHBrCOC,H; 
—> 


COOR COOH OOH 
IV Vv VI 


In this paper there will be described the products obtained from several esters 
(where X =Cl, Br, OCHs;) and benzaldehyde and piperonal. The aryl groups 
are ‘‘tagged’’ so that they can be followed through in reactions used to prove 
structure. Since a chlorinated lactone was first prepared this will be used in 
illustrating the reactions. 

When methyl a-phenyl-8-(p-chloro-benzoyl) propionate and benzaldehyde 
are treated with sodium methylate in absolute methanol, a lactone, C23H,7O;Cl, 
is formed. As the same lactone results when the ethyl ester and ethyl alcohol 
are used, the carbalkoxy group must be involved in the formation of the lactone 
ring. It is possible to write eight possible isomeric structures, of which one 
(XIV) is a beta lactone; the latter is excluded because the product shows 
none of the characteristic properties of this type of cyclic compounds. 


(R2)CesHsCH: OH CsHsCH 7" CsHsCH——-CHCOC,H,X 
| 
(Ri) CsHsC==C——CC,H,X(Rs) CHCH—C—CCHX O=C HC.Hs 
| ta 
o=C—————-O O 
VII VIII IX 
C,H:;CHOH CsHsCHOH CsH,CO 


C,H,CH—C=CC,H,X CoHiC== —— Co.Hs —CH:z—CHC,H,X 
O=C O o=C———-O O= 


Xx XI 
C.Hs,CHOH C,HsCH—O 
CHC —CH=CCH.X c.H.c—cdo CsH,CH:;COCOC,H,X 
sels CHCOGHX 


XIII XIV 
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Any of the y-lactones might equally well be formed; to distinguish between 
them has required a careful study of the reactions. 

Potassium permanganate was rapidly reduced but the oxidation could not 
be controlled and a mixture of benzoic and p-chlorbenzoic acids was formed. 
Chromic acid was the most useful oxidizing agent. In this case benzyl 
p-chlorophenyl diketone (XV) resulted, which at once excluded substances 
having structures XII-XIV. By substituting piperonal for benzaldehyde and 
oxidizing the resulting lactone, piperonyl p-chlorophenyi diketone was obtained, 
showing that it was the aryl group introduced as aldehyde that appeared in 
the diketone. This would exclude a formula like VIII. Further, the substance 
was insensitive to ozone, whereas benzaldehyde would have been easily 
formed from VIII. Finally, the isomeric piperonal lactone (Formula VII, 
R, =3,4-methylenedioxyphenyl, Re=phenyl, R;= -chlorophenyl) has been 
made in this laboratory (21), and on oxidation found to form benzyl p-chloro- 
phenyl diketone (XV). This evidence, taken altogether, excludes all the 
formulas except VII. 

In the Grignard machine (15) the lactone shows one active hydrogen and 
two additions, indicating the presence of an hydroxyl group, but no esters 
were formed in the usual ways, as would be expected of substances X and XI; 
a lactone like IX would probably have no active hydrogen and should add three 
molecules of RMgX. On treatment of the lactone with acetyl chloride a 
chloride was obtained instead of an acetate, a characteristic property of 
tertiary alcohols. Thionyl chloride formed the same halide. This chloride 
(XVI) exhibited many of the characteristic properties of triphenylchlormethane. 
It formed a methyl ether (XVII) with absolute methanol, and an acetate 
(XVIII) with silver acetate. The lactone is thus a lactol. 


C.H;CH: Cl cn OCH; C.H;CH: OCOCHs 


| | 
C.H,C= ——CC,H,Cl C,H,C—=C——_CC,H.Cl a MieaIR das 
0=C———O o=C O o=C 
XVI XVII XVIII 


On alkaline hydrolysis all of these gave a soluble alkali metal salt that on 
acidification regenerated the lactol. 

The latter is insoluble in water and aqueous sodium bicarbonate; it dissolves 
extremely slowly in hot sodium carbonate solution, probably because of the 
sodium hydroxide formed by hydrolysis of the latter. The soluble sodium salt 
is readily converted into a silver salt in the usual way, and the latter gives a 
methyl ester (XIX) when boiled with methy! iodide. The lactol is re-formed 
by alkaline hydrolysis of the ester and acidification. The only evidence for 
the presence of the tautomeric open chain form (XX) is the formation of an 
oxime (XXI). 

C.HsCH: C.H,CH: C.H;CH: 
C.HC=C—COGH.CI C.H.C=—C—COCHCI C.H,C—=¢— CGH! 


COOCHs COOH COOH NOH 
XIX XX XXI 
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A detailed study of certain lactols has been made in connection with keto- 
lactol tautomerism (3, 6, 17, 22, 23, 24) as well as a few instances in suitably 
constituted substances having carboxyl and aldehyde groups in the required 
positions in the molecule (4, 9, 10; 19, 20). In nearly all the cases, typical 
reactions indicating both open chain and cyclic structure were observed, but 
with few exceptions, derivatives of the cyclic forms could only be obtained by 
drastic treatment (e.g., with acetyl chloride or acetic anhydride). A stable 
lactol was formed only with substances having a very highly branched chain 
(XXII) (23). 

OH 
(clt),C-—C—COON 
Oo 


(CH 3)2C—CO 
XXII 


C,H, —" 
C.H;C H=c—¢ HCOC,H; cH—t—CocH, 


OOH COOH 
XXIII XXIV 


The most probable reason for the relative ease with which the lactols described 
in this paper are formed, simply acidification of the alkali salt, is doubtless the 
presence of the very highly branched chain, since very closely related sub- 
stances (XXIII, XXIV) have stabie open chain forms, and only give derivatives 
of a cyclic structure on treatment with acetic anhydride (5, 7). 

The mechanism of the reaction forming the lactol is probably as follows: 
the ester first adds to the aldehyde in the ordinary way, followed by elimination 
of water. The point at which the H and double bond shift is uncertain but 
since Thiele (25) has shown that @, y -esters of this type (X XV) isomerize in 
alkaline media to give a, @-esters the change is exactly what would be expected. 
The water eliminated in the first step in the alkaline solution hydrolyzes the 
ester to an alkali salt; on acidifying the solution the free acid at once cyclicizes 
to the lactol, the tendency to ring formation being greatly increased by the 
highly branched chain. 

CsH;CHO + CsHsCH C.H;sCH; J XX 
C.s,Hs;CHCH,COC,H,CI —> C,H;C H——coGH.Cl C.HC=C—COGH.CI | 


OOCH,; COOCH; OONa VII 
XXV 


The condensation was not brought about by piperidine or diethylamine. 


Experimental 
A. Preparation of the esters. The esters were made as indicated by the outline: 
C,HsCH=CHCOG,H,.X — > C,Hs;CHCH:COC,H,.X 


ry A 7 
+ HCN W i» 


—=—p Se —p- CsHsCHCH,COC,H,X 
COOH “OOR 
‘ D 









CONDENSATIONS OF CERTAIN y-KETONIC ESTERS 609 




















(1) The unsaturated ketones (A). Benzalacetophenone was prepared according 
to the directions given in Organic Syntheses (13). Benzal p-chloro-, p-bromo-, 
and p-methoxyacetophenone were made by the following modification: to a 
mixture of 173 gm. of benzaldehyde, 183 gm. of p-chloroacetophenone, and 
360 cc. of alcohol was added 78 cc. of 10% sodium hydroxide solution, and the 
whole cooled under the tap. The mixture became semi-solid in a few seconds. 
The whole was shaken frequently for 1.5 hr., then the solid filtered and purified. 
The yield was 267 gm. (93%) and the melting point 96° C. The p-bromo 
and p-methoxy homologues melted at 97° and 104° C. respectively. 

(2) The nitriles; addition of HCN (B). a-Pheny]-8-benzoy]-propionitrile was 
prepared by the method given in Organic Syntheses (2). For the homologues 
it was found essential to operate for a longer time and at higher temperatures, 
and to use fresh alcohol in each run; e.g., to 200 gm. of benzal-p-chloro-aceto- 
phenone, 2845 cc. of 95% ethyl alcohol and 87 gm. of glacial acetic acid was 
added 128 gm. of potassium cyanide in 354 cc. of water in 15 min., and the 
whole stirred at 35° C. for eight hours. After standing in the ice box, 193 gm. 
(87%) of a-phenyl-@-(p-chlorobenzoyl)propionitrile, m.p. 122° C., was 
obtained. The p-bromonitrile mixture was stirred for six hours and left in the 
ice box for two days. In preparing the p-methoxynitrile the temperature had 
to be kept at 50-55° C.; both were recrystallized from methyl alcohol. 














TABLE I 
YIELD, PROPERTIES AND ANALYSES OF THE NITRILES 








Nitrile 





Crystal Formula Calcd. Found 
o 


Yield 
% 













p-Bromo- ——— CisH2xONBr Br, 25.5 Br, 25.3 
plates 
p-Methoxy- 65 62 long needles | Ci7HuQ2N | C, 77.0; H, 5.7 | C, 76.8; H, 5.8 

















(3) The esters (D). Some of the esters were obtained by hydrolyzing the 
nitrile to the acid and esterifying the latter (18), and the others directly, by 
saturating absolute methyl alcoholic solutions of the nitriles with hydrogen 
chloride (16), using whichever was found to give the better yield (e.g., with the 
p-chloro derivative, the first method gave a yield of 84% and the second 75%). 






TABLE II 
PROPERTIES AND ANALYSES OF THE ESTERS 





Ester 





Crystal Formula Caled. Found 
form 


| | | || 


Ethyl, p-chloro- 63 fine prisms CigH1703Cl c i. ; Cl, oe 
Methyl, p-bromo- 129 rhombic plates cane Br, 23 0 Br, 22. 
Methyl, p-methoxy- 97 dense prisms CisH gO, C, 72.5; H, 6.0 | C, 72. 3; i 5.9 















h isolated (cf., first method) was not analysed ; by titration with 


The p-methoxyacid (C) thou 
ht of 275 was found (calcd. =284). It formed. small scales from 


standard alkali, a molecular weig 
methyl alcohol. 
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B. Condensation of the esiers with the aldehydes. Asolution of sodium methylate 
prepared from 12.6 gm. of sodium and 212 cc. of absolute methy! alcohol was 
added to a mixture of 106 gm. of benzaldehyde, 212 cc. of absolute methyl 
alcohol, and 160 gm. of methyl a-pheny]-8-(p-chlorobenzoyl)propionitrile, and 
the whole refluxed for two hours. It was then acidified with acetic acid and 
the alcohol and unused aldehyde steam distilled. The residual organic product 
was taken up in ether, the latter dried with calcium chloride, and the solvent 
removed. The solid left was dissolved in the smallest possible amount of 
boiling methyl alcohol and deposited as fine prisms on cooling; the analytical 
sample was recrystallized from an ether-petroleum ether mixture. It is in- 
soluble in petroleum ether, moderately soluble in alcohol and benzene, and 
very soluble in ether and acetone. 

The same lactol resulted when ethyl alcohol and the ethyl ester were sub- 
stituted in the above, and also when the reaction was carried out in alcohol 
that had been distilled from magnesium methylate. The homologues were 
prepared by essentially the same procedure, except that the ethereal extract 
was shaken with a saturated solution of sodium bisulphite when piperonal was 
used. Their properties are collected in Table III. 


TABLE III 
PROPERTIES OF THE LACTOLS 


Found, % 


CHCl 
CoHi70,Cl 


- | 
sos 
a 


2 
7 


¢ 
? 
¢ 
¢ 
? 
¥ 


3 


*From methyl alcohol; > from benzene; ¢ ¢ = C.Hs; 4 pip- = piperonyl ; © from ethyl ester. 


C. General properties of the lactols (VII). The lactols are insoluble in sodium 
bicarbonate solution, but dissolve rapidly in cold, dilute sodium hydroxide, 
and are reprecipitated on acidification. They are recovered unchanged after 
solution in concentrated sulphuric acid. In the Grignard machine (15) they 
reacted with three moles of the reagent and evolved one mole of gas, indicating 
one active hydrogen and addition of two groups. They did not form esters 
(benzoates, p-nitro- and 3,5-dinitrobenzoates) nor a phenylurethane. 
Oxidation. Potassium permanganate in acetone solution oxidized the lactols 
completely to the corresponding benzoic acids. Chromic acid in acetic acid 
did not react as rapidly; a 1,2-diketone was formed and separated from. the 
unoxidized substance by crystallization. Some of the diketones were pre- 
viously known. They were converted into quinoxalines by boiling with 
o-phenylenediamine in alcohol. The lactols (D, F) containing a piperonal 





CONDENSATIONS OF CERTAIN y-KETONIC ESTERS 611 


residue (Re) on oxidation gave diketones that were so sensitive they could not 
be isolated. Even the o-phenylenediamine used in an attempt to get a 
quinoxaline was sufficiently alkaline to destroy them. In all instances, how- 
ever, a deep violet-brown color was given with ferric chloride. A detailed 
description is given of the p-chlorlactol only. 

In a small three-neckéd flask provided with a stirrer, thermometer, and 
dropping funnel, and surrounded by a cooling bath, were placed 20 gm. of the 
p-chlorlactol (VII A) and 150 cc. of glacial acetic acid, and a solution of 6 gm. 
of chromic acid in 15 cc. of acetic acid slowly admitted from the funnel, keeping 
the temperature below 30° C. After an hour the green solution was poured 
into water, extracted with ether and the extract well washed with water and 
dilute sodium carbonate solution. From the latter, on acidification, 3 gm. of a 
mixture of benzoic and p-chlorobenzoic acids were isolated, separated, and 
identified. ‘The ethereal solution, on evaporation, left a residue of 10.6 gm.; 
it was taken up in methyl alcohol, filtered from a small amount of insoluble 
material* and recrystallized to a constant melting point of 103° C. A half- 
gram of the substance and an equal weight of o-phenylenediamine in 15 cc. of 
methyl] alcohol was refluxed for five minutes, and precipitated by addition of 
water. On recrystallization, the quinoxaline melted at 132°C. These proper- 
ties agree with those of benzyl p-chloropheny! diketone (XV) as described by 
Jérlander (11). In asimilar manner the other lactols were oxidized by chromic 
acid to yellow a-diketones. The properties are summarized in Table IV. 


TABLE IV 
PROPERTIES OF THE DIKETONES AND QUINOXALINES 


Diketone Quinoxaline 


| 
Found -D. Calcd. | Found 
% E Formula % % 


plates CuHnOrCl |C,69.6; H,4.2/C,69.5; H,4.2 Ref. 11 

- | needles | CisHuO.Cl Cl,11.7 C1,11.6 CuaHuO:N:Cl | Cl, 9.5/Cl, 9.8 
plates | CiHnOsBr CuHuNeBr | Br, 21.3 | Br, 21.4 
Ref. 11 Ref. 11 


Permanganate oxidation. Several oxidations were carried out, using 5 gm. 
of the lactol in 150 cc. of pure acetone, varying the amounts of permanganate 
and the temperature, but in every instance the result was a mixture of benzoic 
and p-chlorobenzoic acids (4.2 gm. in a typical case) which was separated into 
its components and identified by melting point and mixed melting points. 

D. Properties of the chloride (XVI). The chloride was formed from the lactol 
and thionyl or acetyl chloride equally well. When 15 gm. of the p-chlorlactol 
was dissolved in an excess of the halide, warmed gently, and allowed to stand, 
13.7 gm. of the solid chloride remained after the solvent was removed. After 
several recrystallizations from ether it formed glistening white needles, 

*This solid, m.p. 202° C.; contained chlorine; only enough was obtained for one analysis. 


Found: C, 72.3; H, 4.2%—corresponding to CnHi:0;Cl. Since it was never isolated from other 
oxidations it was not further studied. 
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m.p. 137°C. It is sparingly soluble in cold ether, insoluble in petroleum ether, 
and readily soluble in boiling ether, benzene, and acetone. Analysis: Calcd. 
for CosHO2Cl.; Cl, 18.0%. Found: Cl, 17.8%. In a similar manner the 
p-bromlactol gave a chloride which crystallized in prisms, m.p. 132° C. 
Analysis: Caled. for Ca:H.O,CIBr: Cl, 8.1; Br, 18.2%. Found: Cl, 8.0; 
Br, 17.9%. It was impossible to get crystalline chlorides from the other 
lactols. 

The aceiaie (X VIII). A mixture of 1.2 gm. of the p-chlor-chloride and 1 gm. 
of silver acetate in 15 cc. of absolute ether was refluxed an hour, filtered, and 
the ether allowed to evaporate. The residual solid was recrystallized from 
n-butyl alcohol; it formed shining, white, rectangular prisms, m.p. 157° C., 
insoluble in methyl and ethyl alcohols. Analysis: Caled. for C2sH:Q,Cl: 
Cl, 8.5%. Found: Cl, 8.6%. The same acetate was also formed by boiling 
for 3 min. 3 gm. of the lactol in 10 cc. of acetic anhydride containing a trace 
of sulphuric acid, pouring into ice water and extracting with ether; the yield 
was quantitative. A mixed melting point with the acetate above was not 
depressed. 

Hydrolysis. A mixture of 0.5 gm. of the acetate, 50 cc. of methyl alcohol, 
and 9 cc. of concentrated ammonium hydroxide was left over night at 30° C. 
After neutralizing by addition of dilute acetic acid, crystals of the lactol 
separated. These showed no depression of the melting point when mixed with 
the lactol. 

The methyl ether (XVII). A solution of 11 gm. of the chloride in 75 cc. of 
methyl alcohol was refluxed for two hours. On cooling, the methyl ether 
separated quantitatively, in dense white prisms, and was recrystallized from 
methyl alcohol, in which it is very soluble hot. Analysis: Calcd. for CosHO3Cl: 
CH;0, 7.9%. Found: CH,O, 7.6%. 

Hydrolysis. On refluxing for 15 min. 1 gm. of the ether with 15 cc. of 
10% methyl alcoholic potash and allowing to stand, the potassium salt of the 
acid (XX) separated, and was analyzed without attempting purification. 
Analysis: Caled. for CoxsH»O;CIK: K, 9.4%. Found: K, 8.8%. On adding 
acid to its aqueous solution, the lactol was precipitated and identified by a 
mixed melting point. The ether was not affected by ammonium hydroxide. 

The methyl ether of the p-bromlactol was prepared in a similar manner from 
the corresponding chloride. It separated from methyl alcohol, in which it is 
only moderately soluble hot, as rectangular plates, m.p. 75° C. Analysis: 
Caled. for C4H»O;Br: CH;O, 7.1%. Found: CH;O, 7.2%. It was likewise 
hydrolyzed by alcoholic potash, the ~-bromlactol precipitated by addition of 
acid, and its identity shown by a mixed melting point. 

E. The methyl ester (XIX). The methyl ester could not be prepared by 
refluxing an alcoholic solution of the lactol and a trace of mineral acid (cf. Ref. 
7) but was easily obtained through the silver salt. 

Aqueous silver nitrate was added to 15 gm. of the lactol dissolved in an 
equivalent amount of sodium hydroxide solution, as long as a precipitate 
formed. The latter was filtered, and washed thoroughly with water, alcohol, 
and ether, and dried in a vacuum desiccator. Analysis: Calcd. for Cos3HiO;ClAg : 
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Ag, 22.3; Cl, 7.4%. Found: Ag, 22.9; Cl, 7.6%. A suspension of 18.2 gm. 
of this silver salt in 75 cc. of absolute ether and 10 cc. of methyl iodide was 
refluxed 0.5 hr., filtered, and the solvent allowed to evaporate. There was 
left 7 gm. of ester, which after several recrystallizations from methy] alcohol 
formed short white needles, m.p. 87° C. It was only sparingly soluble in the 
cold alcohol but dissolved readily on being heated. Analysis: Caled. for 
CoH O3Cl: Cl, 9.1%. Found: Cl, 9.1%. ' 

Hydrolysis. Half a gram of the ester was refluxed with 10 cc. of 10% methyl 
alcoholic potash for 15 min., the solution neutralized with acetic acid and 
poured into water. The precipitated lactol was collected, recrystallized, and 
identified by a mixed melting point. 

F. The oxime (XXI). The oxime of the p-chlorlactol (VII) was formed in 
the usual way in dilute alcoholic acetic acid, and purified by recrystallizing 
from dilute methyl alcohol. It formed short needles, m.p. 160° C., moderately 
soluble in cold alcohol and very soluble in hot alcohol and ether. Analysis: 
Calcd. for CosHyOs;NCI: N, 3.6%. Found: N, 3.5%. 
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THE EFFECT OF AGING ON THE ACTIVITY 
OF BAKER’S YEAST' 


By R. K. LARMouR? AND S. F. BROCKINGTON? 


Abstract 


Yeast stored in cakes on ice showed little evidence of change in activity, as 
expressed in loaf volume, during the first 19 Gaps; thereafter the loaf volume in- 
creased until the yeast was 30 days old; after this it decreased but never became 
as low as with the fresh yeast. The rate of carbon dioxide production was quite 
constant up to 26 days; thereafter it increased somewhat and maintained a higher 
oo more irregular rate until the end of the series when the yeast was 56 days 
old. 


Introduction 


Flour is such a complex system of compounds that the task of assessing its 
value and of making accurate comparisons is very difficult and subject to large 
experimental error. Even in the case of a test such as the viscosity test, in 
which the added ingredients are simple chemical compounds and in which no 
manual technique is involved, slight changes in the conditions produce so 
profound an effect that the error of replication is large. The baking test 
involves manual technique which introduces still greater errors due to the 
personal factor of the operator. This factor has been extensively investigated 
by Geddes, Goulden, Hadley and Bergsteinsson (3) and by Merritt and 
Blish (6). 

If all the ingredients used in bread making were simple chemical compounds 
such as water, salt and sugar, the errors observed could be attributed to 
manipulation, but another biological substance, yeast, must be included. 
Thus, there are two uncontrolled factors, the personal factor and the yeast 
factor, and as no one has yet been able to control the former it is very nearly 
impossible to segregate the variability that should properly be allocated to the 
latter. Whether or not the variability of the yeast contributes significantly 
to the total error of the baking test is a question still to be answered. The 
investigations of Werner and Siedhoff (7), Cook and Malloch (2), and Jorgensen 
(5) show that yeasts of different brands may vary considerably. Herman 
and Hart (4), on the other hand, tested two different brands of baker’s yeast 
and found no appreciable difference in the baking results. Thus, all that can 
be said of different brands is that they may vary. 

The product sold under one brand may also vary. Cook and Malloch (2), 
using eight samples of one brand observed the very large range, 289-416 cc., 
in the amount of carbon dioxide evolved in one hour. Another brand showed 
a range of 193-224 cc. for five samples. It is evident that careful factory 
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control may result in a large reduction of variability but, on account of its 
nature, yeast would be expected to exhibit a certain amount of variation, even 
under the most favorable conditions of production. 

Another source of variability is found in the shipment of the product from 
the plant to the distributor, and thence to the baker. This is important in 
western Canada because much of the 
commercial yeast used is produced on 
the Pacific coast, or in eastern Canada, and 
in either case it is several days old by the 
time the user receives it. Although great 
care is used in packing yeast for shipment 
it seemed possible that, in a country where 
wide extremes of temperature are common, 
some variability in activity might arise from 
this cause. An observation made in this 
laboratory confirmed the belief that the 
age of the yeast has a profound influence 
on its activity. In the course of a study 
of variability of loaf volume there was used 
a cake of yeast that had been in storage on Order of baking —S mn. intervals 
ice for eight weeks. The baking results / ; 
obtained with this and with freshlydelivered one flour: A. with neat ecght wecke 
yeast by one baker using the same flour iso 
are shown in graphs A and B, Fig. 1. 

The loaf volumes at the beginning of the series were lower than the values 
previously obtained with this flour and, as the baking progressed, there was 
a very marked decrease in volume until, at the end, the loaf volumes were 
approximately 75 cc. lower than at the beginning. This is not a usual occur- 
rence although occasionally it has been suggested that there was a tendency 
for volume to fall off slightly toward the end of a run, as in graph C, Fig. 1. 
It was very important that definite information should be obtained on the 
behavior of yeast when kept in suspension at 30° C., because in our routine 
baking, yeast sufficient for 25 loaves is suspended at the beginning of the run 
and kept at 30°C. Furthermore, on the basis of the observation of Cook and 
Malloch (2) that yeast kept at 0° C. does not change appreciably in activity 
for periods up to 10 days, the authors had made a practice of having yeast 
delivered once a week and of storing it in contact with ice. It was necessary 
therefore to find out what changes might take place when yeast of varying age 
was kept suspended at 30° C. 


Loat volume ine ¢ 





Experimental 
Series I 
Twelve half-pound cakes of yeast were procured in one lot from the Winnipeg 
distributing plant and placed in cans on ice. Baking tests were made with a 
uniform sample of well aged, commercially milled, first patent flour, and all 
the bakings were performed by one operator. At each baking three doughs 








Loot volume ince 
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were mixed at five-minute intervals with the freshly suspended yeast. The 
suspension was then stoppered and left on the panning bench until the end of 
the baking period, about 3} hr., and shaken at intervals. At the end of this 
time three more doughs were mixed, using the old suspension. Averages of 
each set of three loaves are shown in Fig. 2. In this series of tests, 50 gm. 
of flour per loaf was used in place of 100 gm. During the first 11 days there 
was little change in the volume of loaves made with fresh suspension; with the 
exception of a slight decrease of about 8 cc. on the thirteenth day the volume 
was the same as on the first day, but thereafter a slight rise occurred and the 
loaf volume never fell below the 13-day value. The loaves made with the old 
suspension never equalled the others. The sharp drop in volume shown on 
the eleventh day seemed to be associated with high temperature of the yeast 
suspension. The container had been pushed back too close to the heaters and 
upon examination was found to be at a temperature of 42°C. The previous 
drop on the fourth and seventh days may have been due to the same cause. 
After the eleventh day the temperature of the suspension was checked frequently 
and kept very close to 30°C. After the twenty-first day, however, the loaves 
made with the old suspension dropped in 
volume and remained about 20 cc. lower than 
the check loaves. Disregarding the two drops 
previously mentioned, it appears that in the 
cases of both the fresh and old suspensions 
the activity tended to increase from the 









@ /resh suspension ‘| 
| 2 te off sespasion | 





, - = - + eleventh to the twentieth day, but the old 

Age of yeastin days suspension gave consistently lower values. 

Fic. 2. Effect of age of yeast on After three weeks’ storage on ice the yeast 
loaf volume, first sertes. gave larger volumes than when it was fresh. 


Not satisfied with the results of this series of 
loaf volumes, mainly on account of the error due to failing to keep the 
suspended yeast at a fairly constant temperature, it was decided to repeat 
the experiment with greater care, and to make other tests of yeast activity. 


TABLE | 
EFFECT OF AGING OF YEAST ON LOAF VOLUME, STANDARD V FLOUR 


Loaf volume, cc. Loaf volume, cc. 











Age of | Yeast Age of Yeast 
yeast, | Fresh | suspension yeast, | Fresh  |suspension 
days jsuspension| standing Mean days suspension | standing Mean 
| of yeast | 3 hours of yeast | 3 hours 
| at 30°C. 
| 
7 | 650 29 687 673 680 
12 | 683 33 | 683 710 697 
14 653 41 670 657 664 
19 652 49 | 670 658 664 
26 663 56 | 685 688 686 
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Series II 

For these experiments 12 half-pound cakes of yeast were procured in one lot. 
The number of replicate loaves was increased to five and the basic formula 
with 100 gm. of flour was used. Instead of storing the suspension on the 
bench it was kept in the proofing cabinet at 30+ .5° C. during the interval 
between the first and second mixings. The average loaf volumes obtained are 
given in Table I and shown graphically in Fig. 3. , 

Usually the loaves made with the old suspension were smaller in loaf volume 
than those made with the fresh suspension, but there were two notable excep- 
tions, namely, at 19 and at 33 days. On those dates the loaves made with the 
old suspension were unmistakably 
larger than those made with fresh 
suspension. In the last test, made 


7 








|@ fresh suspension | | 
© Jhr olf suspensien | 





. | Mean 
with yeast 56 days old, there was ; ‘ ee — ae 
no significant difference between 3 : ale . Ne ee” 
the two sets of loaves. These data 3. ° le one 
are inconclusive with respect to occ ® 
the effect of keeping yeast in sus- aS ay | 
pension at 30° C. The average 7 "“Epiduliatlns’ = * 


values were 667 and 663 — for Fic. 3. Effect of age of yeast on loaf volume, 
the fresh and old suspensions _ second series. 


respectively. It may be stated 
only tentatively that in most cases, the old suspension gave slishtly lower 
loaf volume than the fresh suspension. 

In order to compare the trend of loaf volume with age of the yeast, the aver- 
ages of all loaves baked in each test are shown as a curve in Fig. 3, each point 
representing the average volume of 10 loaves. It is unfortunate that there 
are not more data available for the period between 33 and 56 days, as the drop 
and final rise are rather peculiar and may be due partly to experimental error. 
Between 7 and 30 days the loaf volume increased with increasing age of the 
yeast. This improvement in volume was accompanied by improved bread 
characteristics. Interpreting the data broadly, it may be said that after the 
maximum at 30 to 33 days there was a gradual decrease in loaf volume, but it 
should be noted that in no case did the loaf volume drop below the volume 
observed at 7 days. The baking results with yeast varying in age from 7 to 15 
days showed no significant variation and it is concluded, therefore, that for 
experimental baking yeast may be stored on ice for two weeks without under- 
going any change detectable in the baking results. 


Gas Production of Yeast of Various Ages 
Concurrently with the baking tests, measurements were made of the gas 
produced in small doughs. These measurements were made by the method 
described by Bailey and Johnson (1), with certain minor modifications. The 
doughs were made up as follows: flour, 100 gm.; yeast, 3 gm; sugar, 2.5 gm.; 
salt, 1 gm.; water, 63 gm., and when.taken from the mixer were quartered and 
placed in 150-cc. beakers. Two of these were placed in Mason jars containing 
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20% potassium hydroxide solution and two in jars containing 20% sodium 
chloride solution. As rapidly as possible the jars were immersed in a water 
bath at 30° C., and connected to the burettes. Readings were taken at 10-min. 


TABLE II 


EFFECT OF AGING YEAST ON GAS PRODUCTION IN DOUGH. 
RATE OF GAS PRODUCTION IN CC. PER 10 MIN. INTERVAL 


Age of yeast in days 
ime 14 | 19 | 26 | 29 | 33 | 41 | 49 | 56 | Average 
Cc. of CO: produced in 10-min. intervals. 


| 
| 


2 

4 

6 

9 

9 
11 
10 
11 
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11 
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12 

9 
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intervals. This was carried out with freshly suspended yeast and with yeast 
that had been kept in suspension for 33 hr. at 30°C. As there was no significant 
difference in gas production of these two, only the data obtained with doughs 
made from the freshly prepared suspensions have been reported in Table II. 

In dealing with data concerning the production of carbon dioxide, it is 
customary to make graphs of total carbon dioxide produced, volume attained 
and carbon dioxide lost, plotted against time. In this experiment all these 
data were collected but, as we were interested in comparing samples of yeast 
rather than different flours, it was considered that it would be better to calculate 
rates of gas production. Accordingly, in Table II the data are presented as 
the volume of carbon dioxide pro- 
duced in each 10-min. period, and 
thus by inspection the variation 
in rate of gas production may be 
determined. To facilitate com- 
parison of the rates of gas produc- 
tion by yeasts of various ages, the 
average rate per 20 min. has been 
plotted against time in Fig. 4. These 
curves, with the exception of that 
at 56 days, are very similar in 
shape and indicate that the rate of 
gas production in doughs was little 
affected by aging for a period of 49 
days. In all cases the rate increased 
rapidly at first, attaining a maxi- 
mum of 15-17 cc. per 10 min. at the 
end of 70-90 min. Thereafter, the : 
rate of gassing decreased somewhat . Time in min 
rapidly = first and then more Fic. 4. Rate of production of carbon 
slowly until at the end of three dioxide by yeast of varying age. 
hours it had fallen to 10-11.5 cc. 
per 10 min. The curve for 56 days differed considerably from the others in that 
the initial increase in rate was slower and did not reach such a high maximum, 
but at the end of the third hour this yeast was producing carbon dioxide at a 
rate higher than shown by the 7-day-old yeast. 

Probably the most crucial period in the dough is the proofing time, 
which in the authors’ method is between 180 and 235 min. The rates of 
carbon dioxide production at the beginning and end of this period are shown in 
Table III. As the yeast became older there appeared to be an increase in rate 
of gas production during the period corresponding to the proof period. The 
average rate for the 60-min. interval, 180-240 min., was quite constant for the 
yeast samples ranging in age from 7 to 26 days; thereafter, in the 29-day-old 
sample there followed a rather abrupt increase. The last five samples, although 
showing considerable fluctuation in rate of carbon dioxide production during 


Fiate of gas production in ¢.< per 10 min 
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this period on the whole produced gas at a greater rate than the first five of 
this series. The maximum rate for the series was obtained with the 49-day 
sample. 


TABLE III 


RATES OF CARBON DIOXIDE PRODUCTION AT THE TIMES CORRESPONDING 
TO THE BEGINNING AND END OF THE PROOFING PERIOD 


; 
| Rate of carbon | Rate of carbon 
Yeast dioxide production in Yeast | dioxide production in 
sample, cc. per 10 min. | sample, | cc. per 10 min. 
i (from Table II) | age in (from Table IT) | Average 


ays ea 
| | 180 Min. | 235 Min. 





The total amounts of carbon dioxide produced by these samples in periods 
of 3,4. and 5 hr. are shown in Table IV. The values were fairly constant for all 
except the yeasts that were 29 and 56 days old, the former being higher and 
the latter lower than the others. The amount of carbon dioxide produced by 
the 49-day-old yeast was not significantly different from the amount produced 
by the yeast that was seven days old. 


TABLE IV 


TOTAL AMOUNT OF CARBON DIOXIDE PRODUCED BY THE VARIOUS 
YEAST SAMPLES IN PERIODS OF THREE, FOUR AND FIVE HOURS 








Period of Age of yeast, days 
gas 5. 
collection, 





Total amount of carbon dioxide produced, cc. 


309 
$70 | 


z 
238 | 248 | 246 | 223 2s7 | 230 | 226 | 240 191 
| 


319 | 318 | 291 | 333 | 314 | 299 | 309 | 266 


81 | 383 | 356 | 404 | 387 368 382 336 
| 





There is little correspondence between the rate of gas production or total 
amount of carbon dioxide produced, and the loaf volume. The loaf volumes 
increased from 7 to 33 days and thereafter fell off. The rate of gas production 
was fairly constant up to 26 days and then increased and maintained a some- 
what irregular higher level for the remainder of the samples. The total gas 
production was remarkably constant for all samples except those that were 
26, 29 and 56 days old. The mean of the 26 and 29-day samples was in agree- 
ment with the others, but the 56-day sample was decidedly lower. 
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Conclusion 


In conclusion it may be stated that, between 7 and 19 days, yeast stored on 
ice undergoes no change in rate or amount of carbon dioxide production, 
neither is there any significant difference in baking results. Thereafter, the 
loaf volume increases until the yeast reaches an age of 30 or 33 days, after which 
it decreases but the value never becomes as low as the value obtained with 
7-day-old yeast. 

The rate of carbon dioxide production in doughs is quite constant up to an 
age of 26 days but thereafter it increases somewhat and maintains a higher, 
though more irregular, level until an age of 56 days is attained by the yeast. 
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STUDIES IN THE VARIABILITY OF TUBERCLE BACILLI 


Vv. ACID AGGLUTINATION AND ELECTROPHORETIC POTENTIAL 
IN MYCOB. LEPRAE' 


By G. B. REED? anp B. G. GARDINER? 


Abstract 


Previous work has indicated that various species of acid-fast bacteria including 
the tubercle bacilli may be separated into S and R types on the basis of colony 
structure and virulence; some results suggest that the two types differ in respect 
to the surface potential charge on the individual organisms. In this paper it is 
shown that S and R types of Mycob. leprae suspended in distilled water show a 
difference in electrophoretic potential of approximately five times the probable 
error of the determinations. Suspensions of S organisms are shown to have an 
isoelectric point of pH 1.2 compared with a pH of 2.2 in the case of a suspension 
of R organisms. Although acid agglutination of the S and R suspensions was 
found to occur at widely different pH levels, the agglutination occurred at approxi- 
wey, oe same electrophoretic potential for both types, namely, at about 18.2 
mullivoits. 


Introduction 


In an early study of bacterial variation by DeKruif (1) two types of the 
bacillus of rabbit septicemia were differentiated by acid agglutination. Northrop 
(13) and Falk (3, 4, 5) recently reviewed the rather extensive literature, which 
has largely developed since that paper, dealing with variation from the point 
of view of acid agglutination, and particularly with the correlation between 
electrophoretic potential and virulence of strains within a species. It has 
generally been observed that cultures of a particular species showing the highest 
electrophoretic potential show also the greatest virulence. In some instances 
these two characteristics constitute the only observable differences, as in the 
case of the diphtheria bacillus, although this point has been questioned parti- 
cularly by Jones (7). In other cases such as the Pneumococcus there appears 
to be a correlation between potential, serological types, and virulence for white 
mice. Where clear-cut S and R forms have been differentiated it has been 
observed in several instances that the S types, based on colony structure, 
show the higher virulence and the higher electrophoretic potential whereas 
the R forms are avirulent or of low virulence, and exhibit a lower potential. 

Among the acid-fast species it was shown recently by Kahn and Schwarz- 
kopf (8, 9), that S forms of tubercle bacilli show a higher electrophoretic 
potential than the R types. In two earlier papers in this series, Reed and 
Rice (14, 15), it was shown that the S and R forms from a considerable group 
of cultures of tubercle bacilli differ conspicuously in acid agglutination. These 
two results, apparently measurements of the same set of factors, suggested 
the desirability of a more detailed analysis in the hope of uncovering charac- 


1 Manuscript received March 31, 1932. 
Contribution from the Department of Bacteriology, aoare University, Kingston, with 
Sinancial assistance from the National Research Council of Canada 


* Professor of Bacteriology, Queen's University. 
* Graduate student, Queen's University. 
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teristics of variant types of acid-fast bacteria which would lend themselves to 
quantitative analysis. The results discussed in this paper are concerned 
entirely with S and R forms of Mycob. leprae. 

It has just been shown (14) that Mycob. leprae as carried in stock cultures in 
various bacterial collections consists of S forms, R forms, or mixtures of the 
two, together with a variety of intermediate types. In fluid cultures the S 
forms were shown to be highly unstable and to dissociate readily to R. The 
R forms on the other hand proved relatively stable, yet R to S dissociation 
was demonstrated. In addition to the colony structure it was shown that the 
two types differed in habit of growth in fluid media, in oxidation-reduction 
potential changes during growth, and more particularly in complement fixation 
and agglutination reactions with specific S and R antisera. The S type was 
shown ‘to contain a specific S antigenic substance, a species antigen, and an 
acid-fast group antigen, while the R form was shown to possess the two latter 
but to lack the specific S antigen. 


Acid Agglutination 


Well-established strains of S and R Mycob. leprae grown on either solid 
media such as glycerol-egg, or on fluid media such as Proskauer and Beck’s fluid 
were removed with a minimum of media and ground in a mortar with 
distilled water to form an even suspension of the bacteria. This was readily 
accomplished by first grinding the semidry mass from the culture and then 
adding the water drop by drop until an even paste was produced. The S form 
always formed a stable suspension, while the R form required more grinding 


2 TABLE I 
ACID AGGLUTINATION OF SEVERAL S AND R CULTURES OF Mycob. leprae No. 513. THES 
AND R BRACKETED IN THE FIRST COLUMN REPRESENT CULTURES GROWN AND EXAMINED 
UNDER EXACTLY PARALLEL CONDITIONS. THE + SIGNS INDICATE 
THE DEGREE OF AGGLUTINATION 





pH of media 


Organisms 


2.6/2.8! 3.0 


(S 513 
\R 513 


(S 513 
\R 513 


R 513 


+ ++ +/+ +/+ +/+ + 
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++ 44+ 44+ ++ F4+ 
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{S $13 
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and slower additions of water, but where these precautions were taken R sus- 
pensions could be made which remained stable for several hours. The sus- 
pensions were then washed three times by centrifuging and resuspending with 
grinding. Equal volumes of the suspensions and Clark’s phosphate-phthalate 
buffer solutions ranging from pH 2.6 to 7.0 were mixed in small tubes. These 
were incubated for one hour in a water bath at 37° C. and the agglutination 
results read at once, and after 24 hr. at room temperature. The two readings 
were generally the same, the only difference being in the intensity of the 
agglutination. 

Table I indicates the results obtained with several S and R cultures of 
Mycob. leprae No. 513. These were well-established types which had produced 
only S and R colonies, respectively, for several culture generations on solid 
media. The titration results are from cultures of ages varying from 10 to 30 
days grown on gentian-violet glycerol-egg or in Proskauer and Beck’s synthetic 
fluid. The bracketed pairs in the first column of the table represent S and R 
organisms cultured and examined under exactly parallel conditions. It will be 
observed that the pH at which agglutination occurred was uniform for each 
type. 

Agglutination results with several well-established S and R types from other 
cultures of Mycob. leprae are shown in Table II. 


TABLE II 


ACID AGGLUTINATIONS OF S AND R types OF Mycob. leprae. 
THE + SIGNS REPRESENT THE EXTENT OF THE AGGLUTINATION 









pH of buffer solutions 











Organisms ee eh ei ha Sealy RATE wy, Lael ey 
26 | 2.8| 3.0| 3.2 3.6 4.0] 44 | 5.0 | 6.0 | 7.0 

S type 
S 513 + +i+ +] & -~ -~ - _ - — = 
S 65 + ++ 4+ - ~ ~ _ - - _ - 
S 509 + +, - ~ - ~ = - _ = _ 
S 516 + + + - - ~ “ - > sis a 
517 + + - _ _ - - a os - “ 

R type 
R 513 + ++ +/+ +/+ +t ++ +) - - — = 
R 65 ert 4 SF FIT FI OF] - - - 
R 512 + Ht tit Ht +t tit - - _ ~ 
R 519 + ++ +t ++ +/+ 4/14 +1 - | - = - 
R 521 + Ht + ++ 47 FF FH - - - - 


The values stated in Table II represent the average of several determinations, 
as do those in Table I. It is apparent that the S types from the several 
cultures are similar in this respect, and quite different from the R types which 
are similar to each other. 
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Measurements of Electrophoretic Velocity 


Electrophoretic velocity determinations were made with the Northrop (13) 
type of apparatus. In the earlier experiments suspensions of organisms in 
distilled water were used. These were prepared as described in the preceding 
section. As it has generally been observed in such measurements that washing 
the organisms was essential, this procedure was followed in all the earlier 
experiments. It was found, however, that when the organisms were grown on 
solid ‘media and carefully removed, washing did not influence the results. 
All readings were taken at the lower stationary level, as recommended by 
Falk (4), at a distance of }D+4/3 from the lower inner surface of the cell. 
The cell was calibrated as described by Mudd (11). 

Triplicate cultures of S and R forms of Mycob. leprae No. 513 were grown 
for six to seven days on glycerol-egg media, suspended in distilled water, 
washed, and the electrophoretic velocity determined. At least 10 organisms 
were timed in each preparation examined and from these readings the probable 
error was calculated. The results shown in Table III indicate an average 
velocity in microns per second per volt per centimetre of 4.36 for the S type 
and of 3.11 for the R type. The difference, 1.25, is approximately five times 
the probable error of the individual determinations. ‘ 


TABLE Ill 


ELECTROPHORETIC POTENTIALS OF SUSPENSIONS OF TYPICAL S AND R CULTURES OF Mycob. 
leprae No. 513, SIX TO SEVEN DAYS OLD. THOSE BRACKETED WERE CULTURED AND 
EXAMINED UNDER PARALLEL CONDITIONS 


Type of Impressed Number of Average time in Speed in 
organism voltage readings d »/sec./volt/cm. 


10 
10 


10 
10 


10 
10 


Average S =4.36 +.25, 
Average R =3.11+.19, 
Difference =1.25. 


A second group of cultures of both S and R types of varying ages were 
examined by the same procedure. The results tabulated in Table IV indicate 
that the electrophoretic potential increases gradually up to about three weeks, 
and later shows a decrease, especially in very old cultures. The average values, 
it will be observed, are much higher than those indicated in Table III for 
young cultures, however, approximately the same average difference between 
the S and R types is apparent. 
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Use of the Falk Cell 

In an endeavor to simplify the procedure, the depression slide cell as used 
by Falk (4) was tried out. It proved to be less satisfactory than the Northrop- 
Kunitz type of apparatus especially on account of the polarization which 
caused water currents resulting in rapid changes in the readings. However, 
a number of sets of reasonably satisfactory readings are summarized in Table V. 
It will be observed that they are from 1 to 2 u/sec./volt/cem. higher than those 
secured under comparable conditions with the Northrop-Kunitz apparatus, 
but the relative values are very similar and the conclusion the same; namely, 
that the S variant has a higher electrophoretic velocity and therefore a higher 
negative electrical charge than the R type. 


TABLE IV 


MEASUREMENT OF ELECTROPHORETIC VELOCITY OF S AND R TyPEs OF Mycob. leprae FROM 
CULTURES OF VARIOUS AGES, WASHED AND SUSPENDED IN DISTILLED WATER 


S type R type 


Reading in Velocity in Reading in Velocity in 
u/sec./volt/cm. d p/sec./volt/cm. 
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TABLE V 
A COMPARISON OF THE ELECTROPHORETIC VELOCITIES OF SUSPENSIONS OF FOUR TO SIX DAY 
CULTURES OF R AND S TYPES OF Mycob. leprae 1N DISTILLED WATER, 
USING THE F ALK CELL 


SLE ——————— 
Type of | Volts percm. | Time for 420 | Velocity Number of Velocity in 
organism w in sec. u/sec. readings u/sec./volt/cm. 


S type 
S 
S 
S 
NS) 
Average s 


R type 

R | 

R | 

4 
4 
Average R | 
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Determination of Isoelectric Points 


Although the results indicated in the previous sections showed a definite 
difference between R and S types, it seemed from the work, particularly of 
McCutcheon, Mudd, Strumia and Lucké (10), that more striking results might 
be secured by a determination of isoelectric points. For this purpose a series 
of Clark’s buffer solutions, ranging from pH 1.2 to pH 13 were prepared. 
Freshly washed bacterial suspension (2 cc.) was added to 20 cc. of the desired 


S 
S 
3 
x 
X 


' 
‘ 
' 
' 
‘ 
' 
' 
' 
' 
' 
' 


'& Zone of acid 
agglutination 


Fic.1. Curves eA the eno. potential of S pipe curve) and R (lower curve) 


Mycob. leprae iffer solutions of 1topH 7. The data are from 
Table VII :- Seer m X = 10day,O = 14 = O = 15 day, A = 28 day 


cultures; R type; + = 8 day, @ = 11 day, a = 14 day, @ = 20 day, @ = 41 
day cultures. The box enclosed by dash lines indicates the region of acid aggluti- 


nation, data from Tables I and II. 

buffer mixture and electrophoretic readings were taken at once. More con- 
sistent results were obtained when the apparatus was flushed with a portion 
of the suspension in the buffer mixture just before readings were taken. It 
was soon found that there was no basic isoelectric point; attention was there- 
fore focused on the acid end of the pH range. For this reason buffer solutions 
up to pH 7 only were used in most of the experiments. 

Detailed results from the examination of an eight-day R type culture 
together vith the calculated probable errors are shown in Table VI. It will 
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be observed that the electrophoretic velocity varied from 2.01 s/sec./volt/cm. 
at a pH of 5.9 to zero at about pH 2.2. Below pH 2.2 the sign of the charge 
was reversed,and motion toward the cathode became more rapid with increasing 
acidity. Results from the examination of a series of both S and R cultures 
of Mycob. leprae, in the manner indicated in Table VI, are summarized in 
Table VII, the same results being shown graphically in Fig. 1. The R types 
show very consistently an isoelectric region about pH 2.2; in more acid solu- 
tions the sign of the charge is reversed, becoming positive. The S organisms 
on the other hand exhibit an isoelectric point in the region of pH 1.2. These 
characteristics have shown a remarkable constancy in the many cultures 
examined. , 

In the case of mixtures of R and S organisms suspended in buffers of pH 1.5 
to pH 1.7, organisms may be seen to move in both directions. At these 
acidities the S organisms retain a negative charge and move, though slowly, 
toward the anode, while the R organisms assume a positive charge, and move 
toward the cathode. 


Comparison of Electrophoresis and Acid Agglutination 


Ellis (2) found that the stability of an oil-water emulsion was very closely 
connected with the potential difference between the oil drops and the sur- 
rounding medium. Northrop (12), using B. typhosus found a similar condition, 
4.e., that agglutination occurred whenever the potential difference fell below 


15 mv. provided that the cohesive force was not affected. 

It has been shown, Table I, that the R type of Mycob. leprae agglutinates 
with great constancy at pH 4.0 to pH 4.1, while the S type agglutinates at 
pH 2.6 to pH 2.8. In Fig. I, as already noted, the results of a series of 
potential determinations of S and R types have been plotted with the potential 
in microns per second per volt per centimetre as ordinates and the pH of the 
menstruum as abscissa. The pH at which acid agglutination occurs has also 
been indicated on both the S and the R potential curves. It will be observed that 
though the R agglutinate at pH 4 and the S at pH 2.6 to pH 2.8 these points 
are at approximately the same potential, 1.4 /sec./volt/em. Or applying the 


Lamb-Helmholtz equation P.D. = Awe which, as adopted by Northrop becomes 


P.D.=13 Xp/sec./volt/cem., these results then show a critical P.D. of 18.2 
millivolts, which is in substantial agreement with Northrop’s findings that 
typhoid bacilli agglutinate when the P.D. falls to 15 mv. or lower. 

In other words, while the two types agglutinate in the same potential zone, 
it requires a much more strongly acid condition to bring the S type to this 
potential than in the case of the R type. 


Summary 


1. Sand R types of Mycob. leprae have been shown to be distinguishable by 
acid agglutination. The S type agglutinates in buffer solutions of pH 2.6 to 
pH 2.8, the R type in solutions of pH 4.0. 
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2. Electrophoretic potential determinations have indicated a similar type 
difference. With the organisms suspended in distilled water the potential 
difference between the two types amounted to some five times the probable 
error of the determinations. 

3. Isoelectric point determinations have provided more precise and con- 
sistent differences. The isoelectric point of the S types was found to be at 
pH 1.2 and that of the R at pH 2.2. 

4. Although the acid agglutination of S and R types was found to occur at 
widely different pH levels it was also observed to occur at approximately the 
same electrophoretic potential for both types, namely, at about 18. 2 millivolts. 
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A MATHEMATICAL THEORY OF THE GROWTH OF 
POPULATIONS OF THE FLOUR BEETLE, 
TRIBOLIUM CONFUSUM, DUV.' 


By JouHN STANLEY? 


Abstract 


Biological data relative to the growth of populations of the Confused Flour 
Beetle, Tribolium confusum, Duv. have been examined mathematically, the 
individual insects being treated as moving or stationary particles amenable to 
the formulations of the kinetic theory of gases. 

Under certain simplified conditions, é.¢e., prior to the time of the first hatching 
of eggs, it was found possible to integrate the differential equations, obtaining 
curves showing substantial agreement with the biological data. 

Beyond this point, a function @(t) enters, the form of which has not as yet been 
determined, though further work on this point will be carried out. Therefore, at 
present, only a cursory discussion of the use of the function, etc., is given. 

Such information as can be gained regarding the population gro in the later 
stages, without knowledge of the actual form of @(t), is also given. 


Introduction 


The object of this investigation was to work out, as far as possible, a 
mathematical theory to explain the growth of populations of the Confused 
Flour Beetle, Tribolium confusum, Duv., living in whole wheat flour, the cultures 
to be made up and handled as hereafter described. This theory is then applied 
as far as possible to actual experimental data, to show how the various environ- 
mental and biological factors operate to force the growths of the populations 


investigated along the paths which they are found to follow. 

Owing to the fact that all the biotic constants necessary have not yet been 
evaluated, exact numerical solutions cannot be given. Furthermore, owing 
to the complexity of the differential equations, many of them cannot be 
integrated as yet. Hence, the application of the theory must for the present 
be confined to the use of such information as can be obtained from an explana- 
tion of the population trends. 

In order to handle the great number of variables and parameters, a numerical 
subscript notation has been used, whereby it is possible to determine at once 
to which life-history stage, etc., a given symbol refers. 

The Confused Flour Beetle, Tribolium confusum, Duv. occurs commonly 
as a pest of stored products, grain, flour, cereals, etc., in various parts of the 
world. In the United States of America, examinations of such infested 
materials from the more northern states show an infestation consisting largely 
of this species, while in the more southern states, the species, T. ferrugineum, 


1 Manuscript received April 20, 1932. 

This paper was constructed from a thesis submitted to the Graduate Faculty of the University 
of Minnesota in partial fulfilment of the requirements for the degree of Doctor of Philosophy, the 
work being done parily at the University of Minnesota, and partly at the Experiment Station of the 
Association of Hawatsan Pineapple Canners, Honolulu. 
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onolulu, 





GROWTH OF POPULATIONS OF THE FLOUR BEETLE 633 


Fab., is more likely to be found. The two species are almost identical in 
habits and form, and the theories set forth herein can be applied equally well 
to populations of 7. ferrugineum Fab., with suitable changes in the values of 
the biotic constants. 


General Statement of the Habits and Life History 


The adults of 7. confusum, Duv. are small brown beetles on an average 
3.40 mm. long by 1.02 mm. wide, Table VII. The females, which form 50% 
of the adult population, lay small, white, ovoid eggs, at a rate which depends 
on the temperature, Table II. After a period of from 4 to 40 days, Table I, 


TABLE I* 


BIOTIC CONSTANTS.OF T. confusum aT 75% RELATIVE HUMIDITY, AND AT THE 
TEMPERATURES INDICATED 


Coefficient 


Expt. Time in Mean time Standard | Probable of 


days A and B deviation | error variability 


(mean) 


38.80 
38 .83 38 .82 


* The author is indebted to Dr. R. N. Chapman for the data of Table I. 
Tt These are not the A's and B's of Tables X to XIII. 








TABLE II* 


EGG LAYING RATE OF 7. confusum AT VARIOUS 
TEMPERATURES 


Mean , Stand- 


; o- | rate | ard ; Coefficient of 
Temp. °C. perQ | devia- | a | variability 
day | tion | 

TN, ernie 
| 
17 om om ous | on 
22 1.90) 1.181 | 797 | 62.234% 
27 6.24 | 1.067 | 1.124 | 26.622% 
32 | 10.73 | 2.887 | | 1.948 | 26.912% 
i ! | 


* The writer is indebled t to » Dr. R. N. Chapman for the 
data of Table II. 





TABLE Ill 
_ Dimensions oF THE EGGS oF T. confusum * 


































| 
ssa ee al a 
Length | 3 380 | 64 05 
Width | 32 | 80 “40 02 
Data from Brindley (3) (25 measured) tae 
Length | 62 73 | 64 | .04 








Width .38 .47 





40 | 02 


| 


* With adherent flour. 













TABLE IV 
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_MISCELLANEOUS DATA FOR EGGS (STANLEY) 





depending on the tempera- 
ture, these eggs hatch to min- 
ute whitish larvae which 
wander in the flour and, dur- 
ing the course of their lives, 
moult their skins six times, 
increasing in size and vigor 
with each moult. At the last 
moult they change to the 
non-motile pupae. This series 
of changes takes from 17 to 
61 days, again depending 
upon the temperature, as will 
be seen from Table I. The 
pupae remain motionless ex- 
cept for a slight wiggling 
movement for from 5 to 
18 days, once more depend- 
ing upon the temperature, 
Table I, and then, on the 
splitting of the skin, the new 
adult emerges. The new fe- 
male adultdoes not lay eggs for 
a short time, the pre-oviposi- 
tion period, which curiously 
enough is almost the same as 
the duration of the egg stage. 


During life the adults and larvae of various ages wander about in an appar- 
ently aimless way in the flour, the adults, if female, laying eggs en route. 
Both adults and larvae are cannibalistic, living partly upon flour and partly 
upon such eggs, smaller larvae, and pupae as they can find and eat. There 
is no evidence, however, that they purposely search for this living food. They 
appear to be satisfied with that which is offered to them as a result of their 
movement through the flour. They require almost no care to rear, beyond 
seeing that the flour does not become filled with excreted waste. 









Weight of egg, 10,000 weighed (en masse), 
Moisture in mean of 3 samples of 10, 





__ Percentage of eggs to hatch, 27 "=90, U;=.9 
TABLE V 








Stadium 





Length, days | 2.43 3.63 3.03 


m., 0.000,057,8 = W2,(27°C). 
each, 44.958%. 
Available nutrient material in egg, gm. 0. (000,025 =A, 


Note :— Values computed from Brindley's (3) data for 29.7°C, 







STADIA OF LARVAE OF 7. confusum at 27°C. 


Fifth 





Sixth 


3.27 3.39 6.67 
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TABLE VI 
WEIGHT OF ADULT BEETLES OF 7. confusum FROM BRINDLEY (3), 80 OF EACH SEX WEIGHED 


Mean wt., | Stand. 
gm. dev. 


— 


Sex | Min. wt., | Max. wt., 
gm. gm. 

| 
.00006 
.00006 


Male 00140 00155 00148 
Female .00174 .00188 .00178 


From the above general statement, it will be seen that 7. confusum Duv. 
offers certain distinct advantages as a laboratory animal for population studies. 
It makes no webs or other struc- 
tures in the flour, it is fairly re- TABLE VII 
sistant to handling, can be sifted Dimensions oF ADULT BEETLES OF T. confusum FROM 
from the flour to make popula- BRINDLEY (3), 50 OF EACH SEX MEASURED 
tion counts, and can live and ~_ 
develop over a wide range of Dimension Min., Max., | Mean, | Stand. 

mm. mm. mm. dev. 
temperature. For these reasons ————~— — binsetieendieetemtad 
i 
it has been used to some extent Laniah 3.15 3.83 , a 
as an experimental animal, as Width | 85 | 1.11 Po 
the following survey of the liter- ss | 
ature will show. 


A Survey of the Literature 


We shall confine ourselves in this survey to papers dealing with habits, life 
history, population growth, etc., as purely taxonomic papers are not of para- 
mount importance with regard to the problem in hand. 

The most important paper with regard to the insect is no doubt that of 
Chapman (7), in which he shows that the growth of a population of T. confusum, 
is dependent, at any one temperature, upon the size of the environment (a dish 
of flour), and upon the initial concentration of the beetles, but that the final 
concentration is independent of either the initial concentration or the size of 
the environment. Others have checked this work, (See Allee (1), Park (11) ) 
and have obtained substantially the same results. 

That there may be under certain conditions other limiting factors than 
eating is shown by the paper of Chapman (6), in which he speaks of a pungent 
gas given off by the beetle when irritated, and which causes the production 
of monstrosities if present in relatively small concentrations in the air around 
mature larvae. That this limiting factor is not considered in the following 
theory is due to the fact that monstrosities were seldom if ever observed, as 
the populations were always very carefully handled. 

The life history has been worked out by Chapman (4) and by Brindley (3). 
The insect is also mentioned by Chapman (5). 

Its nutritional requirements with regard to vitamins have been studied by 
Sweetman and Palmer (12). 

Further studies have been made upon its life history by Holdaway (9), with 
regard to the production of intersexes, and finally there is the work of Gause (8) 
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in which he shows how changes in the various environmental factors may be 
correlated with the various equilibrium values assumed by the populations. 

Four other papers worthy of mention, as they represent what is up to the 
present the best mathematical work on the subject of population growth, are 
those of Volterra (13, 14, 15), and Bailey (2). 

The writer feels that some criticism may be levelled against Volterra’s work 
on the grounds that so many assumptions have been made in order to simplify 
the mathematical treatment, that the entities considered can nowhere be found 
in the roster of living organisms. 

Although the population counts on which Figs. 2 to 9, and the theory 
itself are based were made, not by the writer, but by Dr. R. N. Chapman, 
and his assistants, an explanation of the technique used would seem desirable. 

Briefly then, whole-wheat flour was sifted through No. 8 silk bolting cloth 
to obtain an homogeneous fine flour containing sufficient vitamins for growth. 
This flour was then placed in the controlled temperature cabinets for a few 
days so that it might come into temperature and relative humidity equilibria 
with the air in the cabinets, and was then weighed out in lots having a weight 
of 32 gm. at 27° C., and 75% relative humidity (air). Previous to this, a 
number of adults were caged and the egg laying rate of each female carefully 
determined. Since violent shaking or sifting was found to alter the egg laying 
rate, these determinations were made by caging single females, and carefully 
rolling them out of the flour with a minimum of disturbance. The eight 
females for each duplicate lot were then selected to avoid the inclusion of any 
females having egg laying rates differing widely from the mean. 

Eight males and eight of the selected females were then placed in each beaker 
with 32 gm. of flour and placed in the cabinets. 

Counts were in general made every 10 days, the various life-history stages 
being separated by the use of bolting cloth sieves of various meshes. It was 
found that repeated counting at intervals of less than 10 days seriously altered 
the egg laying rate. 

The flour was changed at each count, and on returning the beetles to the 
flour, great care was used to distribute them evenly through the mass. 

Thus small controlled environments were set up in which the beetles lived 
and grew, and, provided care was taken to exclude certain parasitic mites and 
intestinal parasites, no trouble was experienced in obtaining parallel results 
with duplicate cultures. It was this close agreement between duplicate cul- 
tures that led to the following investigation subsequent to Dr. Chapman's 
pointing it out to the writer. 


Part 1. 


The Population Growth Under Conditions Such That 
Environmental and Biotic Resistance are Zero 


It will be apparent that differential equations can be written, descriptive 
of the growth of the population, by the use of parameters descriptive of the 
insects and the environment. 
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In the following discussion the symbols pertaining to the various materials, 
flour and the different life-history stages, will be identified by the use of the 
following numerical subscripts, in accordance with the scheme outlined in the 
introduction. 

Flour, 1; eggs, 2; first instar larvae, 3; second instar larvae, 4; third instar 
larvae, 5; fourth instar larvae, 6; fifth instar larvae, 7; sixth instar larvae, 8; 
pupae, 9; immature adults, 4.e., adults which have not passed the pre-oviposi- 
tion period, 10; mature adults, 11. 

Since however the point which we desire to make in Part I can be proved 
without a consideration of the various instars, etc., there is given below only a 
limited notation, sufficient for the matter in hand. 

Let: No=original number of adults. 

N2=number of eggs at any time, T. 
R=sex ratio, t.e. ratio of females to total. 
€=number of eggs laid per female per unit time. 
To=origin of time, 7.e., time at which population is set up. 
t2(M)=time at which the laying of eggs of the M‘ generation com- 
mences. 
ts(M) =time at which the hatching of eggs of the M‘™ generation com- 
mences. 
I’2=number of days spent in the egg stage. 
S= total time from egg to egg. 

We shall now, under Part I, consider the idealized growth of a population, 
i.e., on the supposition that there is no environmental resistance, no can- 
nibalistic eating, and no variation in the values of the biotic or environmental 
parameters. 

It will be seen, after careful consideration of the matter that, in any genera- 
tion, there are two phases. For example, in the M‘" generation there are the 
following:— (a) Phase M,—from the commencement of laying by the M‘™ 
generation adults to the time of hatching of their eggs, 4.e., during the period 
of time, #(M)<T<4(M). (b) Phase M,—from the commencement of 
hatching of the M*® generation eggs to the commencement of laying of the 
(M+1)™ generation eggs, that is during the period, 4(M)<T<4(M+1). 

The equations descriptive of the growth of the egg population during the 
various phases are as follows. 


Phase dq; 


This phase covers the time interval, #.(1) < T< (1), ¢.e., from the commence- 
ment of laying of first-generation eggs, to the commencement oftheir hatching. 

During this period, 

N,=ReENoT 

Where Np is the original number of adults, not necessarily 16, as in the 
experiments performed in connection with this particular problem. 
Phase 1, 

This phase covers the time interval t;(1) < 7< (2), 7.e., from the commence- 
ment of hatching of first-generation eggs to the commencement of laying of 
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second-generation eggs. During this phase, and during all subsequent time, 
the rate of increase of eggs is diminished by the hatching of first-generation 
eggs. The number of first-generation eggs which have hatched up to a time T 
is the number which were laid up to a time (7 —J>), that is, 
ReNg(T —I>). 

Whence, during Phase 1; 

N2 = RENoT — RENo( T —I2) = RENoI2 = constant. 
Phase 2, 

This phase covers the time interval, /.(2) < T7< (2), ¢.e., from the commence- 
ment of laying of second-generation eggs to the commencement of their 
hatching. During this period, the rate of change of N,2 is increased by the 
production of second-generation eggs which are laid by adults formed from 
first-generation eggs laid during the period 7)>< ¢< T—S, where S is the time 
necessary for transformation of a newly laid egg to a fully mature adult. 


The group of NoR first-generation eggs laid at T= 7+ — will mature and 
lay eggs in turn at T=S+ and will lay NoR?(T—S)e second-generation 
eggs up to a time 7. 

The next group, laid at T=S+ < will in turn produce second-generation 


eggs to the number of N,R*[(7T—S)e—1] up to a time T. Thus the total 
number of second-generation eggs laid subsequent to a time, S, and up to a 
time T, will be, including the last group, 


wana) (rset [P-Spe—t] +[(7-Se—2]+ ee (7-5-7 


a =(T—S)e 
= YL MR|(T-S)e-a). 
a =0 
Whence, during Phase .2,, the total number of eggs at any time 7, during the 
period of the phase is, 
a =(T—S)e 
N:=ReNI3+ LD NoR*[(T—S)e—a]. 


a,=0 
Phase 2; 

This phase covers the time interval, t;(2)< 7<4(3), ¢.e., from the com- 
mencement of hatching of second-generation eggs to the commencement of 
laying of third-generation eggs. During this phase and after, a negative term 
enters due to the hatching of the second-generation eggs. The number which 
hatch up to a time T is equal to the number laid up to a time (7 —T>;), that is, 

a = (T—S—T se 
> NoR?[(T-S—I2)e—a ]. 
a =0 
Whence, at any time, 7, during this phase, 
a = (T—S)e a =(T-S—Te 
Ni=ReNI1+ = = MR (T-S)e—a]—- E NoR?[(T-S—I2)€-—a J. 


a, =0 a,=0 
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Phase 3; 


This phase covers the time interval t.(3)<7<4#(3), ¢.e., from the com- 
mencement of laying of third-generation eggs to the commencement of their 
hatching. Ina manner analogous to that in which the second-generation eggs 
discussed under Phase 2, were formed, these third-generation eggs are formed 
from eggs laid during a time subsequent to T= 25S. 





The second-generation eggs laid at a time 7 =2S+ “ are equal to NoR’, and, 


upon transformation to adults, will lay third-generation eggs to the number of 
N,R®[(T—2S)e], up to a time T. 


The second group of second-generation eggs laid in the interval 
i=2 
2S+ ; <T<2S+ . are, from the calculations of Phase 2;, equal to = R% 


t= 


2 
and lay third-generation eggs to the number of } NoR%|(T—2S)e—1] up 


i=l 
to a time T. 
Whence, by a continuation of the above reasoning, the total number of 
third-generation eggs laid up to a time T is, 
i=1 i=2 
N=  N R4u[(T-25)e-0]+ DY MRi[(T—25)e—1]+ 
t=1 i=1 
i=3 


D NeRi[(T-25)e-2]+ .. 2.2... 
t=! 


i=(T—29)€ 1 
Dicuss Za NoR'i[(T—25)e—(T-2S— —)] 
i=l 
@,=(T—2S)€ Mea: 
whence, N= > DY MR {(T-25)e—an J. 
Oy=1 a =0 


Whence, during the interval, t.(3)< T<4;(3) 


a =(T—S)é a =(T—S—T)) 
M=ReENI2+ =D = NoR*|(T—S)e—a|— = NoR?[(T-S—T2)€—an | 


a, =-0 a, =0 










G@:=(T—2S)€ ima 


eZ 2 NoR*[(T-2S)e—ex:| 


a@,=1 a= 
Phase 3: 

This phase covers the time interval, t;(3)< 7<4#(4), 4.e., from the com- 
mencement of hatching of third-generation eggs to the commencement of 
laying of fourth-generation eggs. By an extension of the reasoning employed 
under Phase 22, the additional negative term is 








@.=(T—S—To)€ a =a 


i > > N.R*[(T—S—I2)€—a |. 


a@,=1 a, =0 


Phase 4, 


This phase covers the time interval, &(4)< 7<4(4) i.e., from the com- 
mencement of hatching of fourth-generation eggs to the commencement of 
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laying of fifth generation eggs. By an extension of the reasoning used under 
Phase 3,, the additional positive term is 


@&=(T—-3S)€ Giea; aa 


z NoR‘[(T-3S)e—a J. 


@=1 Qy=1 Oh = 
Phase 4, 

This phase covers the time interval, 4,(4)< 7<#(5), 4.e., from the com- 
mencement of laying of fifth-generation eggs to the commencement of their 
hatching. Again a negative term enters, which can be shown to be 

@=-(T—S-I)¢€ =a: M=mH 


a ZZ MRi[(T-35-M)e—ar]. 


y= &=-1 He 
Whence, for phase M; where & (M)<T<4(M) 


a =(T—S)€ (T-—S—faeé 
Ni=ReNI+ LY = NoR[(T-S)e—ai]- 2, NR*[(T—S—I:)e—an] 


a,=0 aQ- 


One (T—2S)e aids @=(T-—2S—Fie aay 
+ ZX NR {(T-25)e—an] — N.R*[(T-2S—Ts)e—ax | 
e 0 


a,=1 aun a@,=1 a= 


Myg—1) M(T—(M—1)SJ€ Bi yg—2) = Oa —1) 
oa 


@y—1) 1 


+ > z wan] fr (ue-t)s hema | 


a@y=1 a= 


During Phase Mz, where é;(M)< T<#(M-+1) the additional negative term is 


© Uag=aj ~(T — (64-1) S— Fife eg - 2) = % 44 -1) Oi -) %-84 1 


=1 1 a =1 


Ou -2)™ 


e & “S- ner” [ {r-ae—ns—ripena |. 


ae=-1 a@=0 


* im —1) (M-® 
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By an extension of the above method we may write equations descriptive 
of the growths of populations of any life-history stage. 

The above equations demonstrate two important facts: 

(a) That a population of this type, having distinct generations does not 
increase according to the compound interest law. It could probably be shown 
that such a law of increase does hold when T becomes very great, or if S is 
very small, as for instance in the case of bacteria and protozoa. 

(b) That the population existing at any time is a function of its whole past 
history. Hence damage to a population is, theoretically at least, irreparable 
within finite time. In actual populations, owing to the action of environ- 
mental resistance, damaged populations do tend to return to the form which 
they would have had if undamaged, but never actually reach it. 


Part 2 


Population Growth Under Conditions Such That Environmental 
and Biotic Resistance are Greater Than Zero 


The population growth will now be considered under the assumptions that 
temperature, relative humidity and the weight and chemical constitution of 
the flour are constant, while the various subtractive forces, resulting from 
natural mortality and from the eating of one form by-another, are allowed 
full play. 

These assumptions bring a very complicated system of forces into existence, 
hence the following general statement, supplementing that already given, 
may be of interest. 

The original 16 adults, eight males and eight females, are placed in the 
flour and at once commence to lay eggs. At the same time, as a result of 
their movement through the flour, they again encounter some of their own 
eggs, and may eat them, certain conditions being satisfied. As a result of 
this eating and Jaying, the eggs tend to increase to the point where they are so 
numerous as to be found and eaten as rapidly as laid. Flour is of course eaten 
at the same time, a certain mathematical relationship existing between the 
amounts of egg material and flour consumed per adult per unit time. 

Subsequent to hatching, the new larvae themselves prey upon the eggs, 
consume flour, and are in turn preyed upon by each other and by the adults. 
Young larvae also fall victims to older larvae. 

At the time of the first pupation, the larval population is temporarily 
decreased by transformation to the pupal stage. At the same time, the egg 
population rises owing to the fact that egg-eating larvae are being transformed 
to non-egg-eating pupae. Thus at this point, except under certain special 
conditions to be discussed later (see page 666) the egg and pupal populations 
are increasing, and the larval population is decreasing. 

Upon the commencement of emergence to the adult form, the passive pupae 
are changed over to voracious egg-eating adults, as yet too immature to lay 
eggs. Consequently, egg, larval, and pupal populations tend to decrease, 
while the adult population increases. 
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After a certain time, the pre-oviposition period, the oldest of the new adults 
z.e , the oldest female, begins to lay, and from this point on, the egg population 
increases to come into periodic fluctuating equilibrium with all the other 
forms. This increase in number of eggs is reflected in a later increase in 
larvae, and a still later and smaller increase in pupae. Occasionally there 
may be a later and small increase in adults. The eating of one form by another 
exerts a damping influence on the amplitude of these small increases. Figs. 4 
to 9 will show the general truth of the above statements. 


The problem of the numbers of contacts in unit time, occurring between insects 
moving in the flour 

In order to discuss the manifold forces acting within the population, it is 
necessary to have some means of computing the rate at which contacts occur 
between the various individuals moving in the flour, since such contacts may 
result in the eating of some individuals by others. 

If the individual insects be considered merely as moving particles, and 
owing to the almost infinitesimal degree of intelligence possessed, it is legi- 
timate to make such an assumption, the problem may be handled by means of 
formulas from the kinetic theory of gases. 

This theory (10) shows that three types of contacts may occur between 
particles moving in a given space: (a) contacts between a moving particle and 
a stationary particle, (b) contacts between two particles moving with the same 
speed, and (c) contacts between two moving particles having different speeds. 

Thus, under case (a), where V is the volume of the space through which the 
particles are distributed, v is the number of stationary particles, o, their mean 
diameter, 0, the diameter of the moving particle, and » its mean speed, the 
number of contacts in unit time between the moving particle and the » station- 
ary particles is, 
oi +o \? 
ver (242) 

Na= a (1) 


Under case (b) where both particles move with the same mean speed yz, 
and have the same diameter o, and where the total number of particles per 
cubic unit is vy, the numbers of contacts per unit time between any one moving 
particle and the remaining (v—1) particles is, 

+ ux(v—1)o* 
Bye = some (2) 
V— 3 (v—1)e 


Under case (c), where the two types of particles have mean speeds s; and fz, 
mean diameters of 0; and o2, and where the number of one type is vy, the number 
of contacts in unit time between a single particle of the other type and the v 
particles is, 

1 


in 641 me 


aa vo Es) 
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The value of r= ia [ = (er 1s)*-+ (eto)? | depends on whether 2 > #2 or 
fo > pr for if 41>. then, 


Sun? + ue? . 
3m 


and if #2> then, 
3 yo? + . 
3 2? 


This difficulty may be overcome by writing 
_— Guus [ —(|m—ml)? + (+00) | = F(u, me), 


where the vertical lines around “4—/2 have their usual significance, meaning 
the absobute value of 1: — fe. 


T= 


In the case of contacts between individuals of the various life-history stages 
of T. confusum, it will be apparent that larva-egg, larva-pupa, adult-egg and 
adult-pupa contacts come under case (a); that larva-larva (where both larvae 
are of the same age) and adult-adult contacts come under case (b); and that 
larva-larva (where the two larvae are of different ages), and adult-larva contacts 
come under case (c). 


2 
Since, moreover, the volume V—$zv (257) is merely the volume, G, 
of the flour, the above formulas reduce to 


(2 77) 
ME 2 


G 


Na= 


4 ux (v-1)o? 
SS ee ae eee 


G 


F(mi, wa) ev (“77/ 
N.= a nen (7) 

It is necessary at this point to discuss more fully the exact significance of 
0; and o2 in the case of insects. It is not necessary that an insect come into 
actual physical contact with the main body wall of another insect to be aware 
of its presence. It may for example touch against a minute projecting bristle. 
As, then, a ‘contact’ will be defined as any approach of two entities of the 
‘ population within such a distance of each other that mutual recognition of 
each other’s presence occurs (or recognition by one party occurs, if the other 
be an egg or a pupa), it will be necessary to define what may be called “radii 
of perception”’. 

It has been suggested to the writer that since the insects under discussion 
are assumed to have no intelligence, the term radius of perception is somewhat 
unfortunate. It is felt, however, that such is not the case, since by 
“intelligence”’ is meant, not the mere possession of powers of perception through 
the senses, but rather the possession of mental powers sufficient to make 
conscious and more or less elaborate decisions under definite conditions. 
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The radius of perception of a living entity for an infinitesimally small point 
is defined, then, as equal to the radius of a sphere whose volume V is 


8 Yi 2 
V= [ I [ wx, f, Z)dx dy dz, (4) 
a Ye JZ 


where J(X, Y, Z) is a surface surrounding the entity in space, beyond which 
its powers of perception are zero. 

However, neither party in an insect contact such as have been discussed 
above has an infinitesimal diameter. Hence, where, for example, the radius 
of perception of an adult for an infinitesimal point is 7,2 and the radius of 
perception of a larva of some certain instar such that the correct subscript 
may in a general way be written as 7, for the point is, 7;,, the radius of percep- 
tion of either the larva for the adult, or vice versa is 


Ty Ptthiy 
fs? ——— 


In the case of contacts with eggs, the latter of course are incapable of per- 
ceiving anything so that, strictly, r,;=0. It is obvious however that the 
quantity to be used in such a case in place of 7, , is the mean radius of the egg. 

Then where M,,,, My, and M,,; are the numbers of contacts in unit 
time between adults and eggs, between generalized larvae of the same age, 
and between adults and generalized larvae respectively, and where p,; and pz, 
are the mean speeds of adults and generalized larvae, 


aMiy2 _ Nu Nouss2ri1.2 
a G 


dMir _ 4N(Np—1)un rh 
at 3G 


aMy 1 = Nyy N° F(uss, wn) ®rin.2 
dt G 


Amounts of Flour Encountered by Adults or Larvae 


Since adults and larvae eat flour, as well as eggs, other larvae and pupae, 
it is necessary to develop functions descriptive of the rate at which flour is 
encountered by moving adults or larvae. Theoretically this function could be 
developed in a manner similar to the above-mentioned formulas for contacts 
between adults and eggs, etc., but a much simpler and neater derivation can 
be given by reason of the more or less continuous nature of the flour medium. 

It might be felt that inasmuch as an adult or larva is continually surrounded 
on all sides by flour, and as, by reason of its homogeneous nature, one portion 
of flour is indistinguishable from any other portion, the amount of flour 
brought to the attention of a moving insect in unit time should be independent 
of the mean speed of the insect. (It should be noted that in speaking of the 
senses of insects, no anthropomorphic concept is intended; the words are used 
simply because more suitable expressions are not available.) 
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The writer feels however that such is not the case, for he cannot help but 
think that (to use unfortunately an anthropomorphic analogy) a steady stream 
of material passing across the field of perception gives an impression of volume 
proportional to the distance travelled by the stream in unit time. 

Consider an adult moving in a random way in a space, G, filled with flour. 
Clearly the adult does not bore out an exactly circular tunnel but, for ease of 
consideration, the tunnel will be considered as circular, and of such a radius 


p that zp?= [ xdy where the expression on the right is the line integral taken 
c 


in the positive direction around the boundary curve of a right section of the 
tunnel. 

It is further assumed that the anterior surface of an adult cephalad of a 
certain plane of right section Q is defined by a function ¢ (x, y, 2); that only 
the flour in a layer of thickness D over the surface ¢ (x, y, 2) is perceptible to 
the adult as food, and that when, due to forward movement of the insect, a 
particle of flour passes to a position posterior to Q, it ceases to be perceptible 
as food. It is also assumed, for the moment, that flour is a perfect fluid, i.e., 
frictionless and incompressible. 

If now, the surface ¢ (x, y, z) be moved a distance dx along the longitudinal 
axis of the insect, a certain frustum of thickness dx will pass beyond Q and 
be lost to perception. 

With the exception of infinitesimals of higher order, the volume of this 
frustum is: 

dv = x(9+D+p)(o+D—p)dx 
= x(20D+D")dz. 

Whence, on advancing the surface ¢ (x, y, z) at a rate #, there passes across 
the field of perception of an adult, in unit time, 

Hu 


V= | «(20D +D")dx = x(20D+D*) yu), 


which is the volume of a cylindrical shell of thickness D surrounding the bore 
of the tunnel cut out in the flour. 
Then, where, 
M,,., =volume of flour encountered by a mature adult up to a time T. 
M_, =a similar function for the generalized larva. 
And where, for the sake of uniformity, we write: 
ria =20uDutDi ri. =2p~Di+Di 
Calling r,,;, and r,, the radii of perception for flour, we have: 
dMi1 1 2 aM, 2 
SS = 877.1 BL: 
a 
In a practical case, however, flour is not in any sense a perfect fluid, but this 
makes no difference as far as the speeds 1, and yu, are concerned, as these are 
the speeds maintained in spite of the resistance of the flour. The fact of 
compressibility must however be taken into account, and this may be done by 
means of an arbitrary factor of compression, to be determined by experiment 
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for each type of flour and for each life-history stage. These factors are then 
of the form, K3, Ky, Ks, Ke, Kz, Ks, Kio, Ku, according to the established 
notation for subscripts. 
Including the above factors, where K, is the subscript for the generalized 
larva, 
4aM444 aM, 


1 =Kyer ; 1 =K 2 : 12 
4 171.114 a LITLABL (12) 


A Transformation of the Above Equations 


So far the equations descriptive of rates of contact have been developed in 
terms of individual eggs, larvae, etc., but in order to determine the rates at 
which these various food materials are consumed it will be necessary to trans- 
form the above equations into others descriptive of the rates at which assimil- 
able nutritive materials are encountered. 

It would be difficult or even impossible to determine the exact requirements 
of each life-history stage in terms of the various fats, proteins, vitamins, etc., 
and if the various stages were at certain times dependent entirely upon one food 
and at other times entirely upon another, it would be necessary to have exact 
knowledge of the materials obtainable from each food source. However, 
as this is not the case, it is possible, without grave error, to make the assumption 
that each food material supplies all of the substances necessary to maintain 
the various entities feeding on it. It may be the case that eggs, for example, 
are lacking in some particular chemical substance if they are to function as the 
sole food, but this need not cause concern, as the necessary compound is 
obtainable from some source as evidenced by the vigorous growth of the feeders. 
That is to say, it does not appear that, throughout all the available foods, any 
compound is scarce to the point of detriment to the population. 

Suppose that A;, As, As; etc., are the assimilable percentages of nutrient 
material in the various life-history stages, W2, W3,etc., the weights of individuals 
of the various stages, and W, the density of uncompressed flour. We obtain, 
then, as equations descriptive of the rates at which assimilable nutrient 
materials are encountered, the following : 


d¥ i342 A,W2Niy Nous 711.2 
at G 


€¥inc _ ALWiNuNz* F(us1mz) 81.2 
dt G 
aY;2 _ A,W,N, Nou ery. 


dt G 


_ AW Kin eri1aG 


G = AWW Nu Kiswu tris 


© AWNLK pati sG 


G =A LWLN LK perry. 
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It is now possible to set up the differential equations for the early stages of 
the population growth, except for one point discussed below. 

It will be apparent that (considering the insect from a purely mechanistic 
viewpoint) there is a certain definite quantity of nutrient material which must 
be assimilated in unit time to carry on the processes of life during that unit 
time. (The question as to what would occur in vitamin-deficient foods need 
not enter here as the use of whole-wheat flour obviates'the possibility of such 
deficiencies.) 

In actual practice it is not certain that the various life-history stages feed 
each one at a constant rate throughout the day. However, as they were kept 
in darkness, and the flour environment is at all timés invariant, and there is at 
all times an abundance of food, it is assumed, lacking evidence to the contrary, 
that assimilable materials are ingested at such an average and approximately 
constant rate that the maintainance ration of the same is obtained for each 
increment of time. It is apparent that such an assumption cannot exactly 
describe what goes on, as, if larvae of an age a happen to be singularly lacking 
in assimilable materials, an adult, while consuming one of them, cannot, 
perhaps, maintain this assumed rate of ingestion. However, this is not a 
grave matter, as the above assumption means that, with the exception of 
negligible quantities, the required amount of material is ingested over a finite 
time, even though the rate of ingestion may be small at some one instant. 

This assumed, approximately constant rate of ingestion will be referred to 
as the “maintainance rate of ingestion”, with the additional assumption that 
the digestive powers of all life-history stages are the same. Symbols of the 
form E;, Ey, etc. will be used to denote it. 








The Population During the Period t,=f.(1)< T<t#;(1) 


During this period, mature adults are laying eggs and at the same time are 
finding and eating them, whence the rate of change of the egg population is 
equal to the difference between the rates of egg production and egg con- 
sumption. 

Let: 

Cis.2 = number of individual eggs consumed by the N,,; mature 
adults up to a time T. 

Xi1.2, Xi1.1=amounts of assimilable material, (by weight) obtained 
from eggs and flour respectively, and consumed by the 
Ni, mature adults up to a time T. 

Pi:.2, Pi:.,=constant coefficients of preference of mature adults for 
eggs and flour respectively, as foods. 


aN 
Then: = = ReNy, —C's3.2 (18) 














From the definition of P,; 2 and Py; ;, 
—_X'n2 Xing 
Pia Y's2 Pass Y' tas 





(19) 
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From the definition of En, 
Mp2 X 1 = EN i- 
From (19) Xigg we ee 2 
Pii2¥'n12 
Substituting in (20) 
Pua Yuna X12 


sas = : = Ey, 
Pu2¥'s2 


y yr 
—EaN Pua u2 


X'n2= = , 
Pay 2¥' spat Piva Y's 


Ey NyPi1.242Wamrr erin aNiNa 
G 


Py 2A.Womn erin aun + PsA Ws Kumari GN 
G G 


X'n2 _ EuPyarin2NiuN2 


C'13.2 = a ee dee oe ee ae ’ 
AgW2 P43 2AaWoris 2No+P iy sArWiKutinaG 


(24) 


HIN, 
cN,+4 
Where H=Ey Pyy2 ira Nu; c=Pyi2 Ae Wa riva3d=Pys Ai Wi Ku tis G. 

Whence 


Cur (25) 


aN, ReN,,cN2,+RENid—-HN, 


(26) 
dl cN,+d 


aN, aN2z+b 
dat cN,+d 


Where a=ReNi,c—H; 6=ReNid. 


Separating the variables and integrating: 
cN,  ad—be 

a + a2 
When t=0, N;=0, 


So that _ad—be 


‘= 


log(aN2+6) =t+a 


log 6, 


a2 
c ad —be aN2+b 
T= @ Nat ae log (+) = F(N:). 


Whence 


Characteristics of the function T = F(N:), (See Fig. 1). 
We may write, 


aN, +h = (ReNy,6+H)N2 +ReENj,d 
ONg+b = NyP43.2711.2|R€AgW2—E|No+ReNysPy1,:41WiKurin.sG, 
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from which it is evident that a is less than zero in any real population, since 
ReA.W2, the amount of assimilable nutrient material expended per female 
per day in the form of 

eggs, must be less than ¢ 





Ey, the total amount 
of assimilable nutrient 
material taken in per 2 L 
female per day, i.¢, & : 
(ReA2W2—En)<0. ‘ i 
Hence, as b, c and d are 
easily seen to be greater inieiaatiiiiligigilaiaiad 
TIME 
than zero, 
Fic. 1. Graph of the function, T = F(N:). 
dN,> < 6 
a EO. a em ae (30) 
ad < = a 


It can further be shown that T=F(N2) may be written as a convergent 
alternating power series of the form, 


No =a, T —apT*+03T?—ayT*+ 2.2... , (31) 
and that a =ReNy; ® 
whence lim. Nz 


TO ae =ReNi, 





Also, since dN, aN,+b lim. dN, 6b — lim. aN; 
os ’ —— of we + = — . 
a: ee Ore Ss . 


It is also evident from an examination of (31) that aan N, =0 and it can 


‘ : ee lim. @N. 
further be shown by direct differentiation that ican = =0 





It is thus apparent that the egg population is equal to zero when T is equal 
to zero, and increases thereafter to approach a value £ which it reaches only 
after an infinite time. (In an actual case, owing to the fact that the function 
T = F(N2) can exist only for integral values of N2, this limit, §, may be reached 
in finite time, the population oscillating around the value &, and between the 
limits of the least integer greater than £, and the greatest integer less than 
g.) 

Furthermore, when 7 =7>, the egg population is increasing exactly at the 
rate at which the N,, mature adults lay eggs, namely, at the rate ReNu. 
This rate decreases thereafter, reaching zero when N; equals £. 

It is interesting to note the meaning of —. It represents the point at which 
the eggs are found and eaten as rapidly as laid. Thus, if eggs are introduced 
into the system by artificial means, as by adding them to the flour by hand, 
and stirring them in, they will be found and eaten more rapidly than the Ni, 
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dN, ° 
adults can lay, so that 37. <0 and the egg population decreases to approach 


the limit — from above. Such conditions as have been described above occur 
only, of course, in the absence of hatching. 

It is also of interest to examine the relationship between ¢ and G, the size 
of the flour mass. Since hatching does not occur, Ny, is a constant. 


We may write, b = ReN,,d = ReN, Py :41Wi Kunis iG = yG, 


whence YG 
= a ’ 
and i. = —— ~<a constant . (32) 
dG a 
Whence, the final egg population (in the absence of hatching), is proportional 
to the size of the environment. Since we may further write, a=0Nu, 


where @ = Py, A2Woris.2(REA2W2—Ej) 

and b=ONy , 

where $= ReAyW,KyrinaG 

Se a Ee. (33) 
ON 4 

which is independent of the value of Nu. 

Chapman (7) has shown that the limiting population reached with T. 
confusum is, by actual experiment, proportional to the size of the environment, 
and independent of the initial number of adults. It is true that at any temper- 
ature above 17°C. hatching occurs, so that the final population is a mixture 


CL Lie eo ie ee ed 
od eee CL 
Ld ene HFS # 
| | aoeroom—o! TE AL OS 

Shenae == ad ash aaa 

LJ Coe | | 
| ttt ttt ttt Te 

hhh het tht teh heetemfimk | 


° 
0 0 20 30 #0 30 60 70 80 90 0 110 120 130 140 150 160 170 180 190.200 210 eee 
TIME IW DAYS 


Fic. 2. Population growth of T. confusum, Duv. at 17°C., mean of (A) and (B) populations. 


Then t= 


ote 
oO 0 20 30 % 30 G0 70 Bo G0 700 7 130 140 150 160 170 1B) 100200 eseues) 
TIME tl DAYS 


Fic. 3. Smoothed curve for population growth of T. confusum, Duv. at 17°C., (See Fig. 1). 
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of all stages, but it can at least be said that theory and experiment are in sub- 
stantial agreement up to 17°C. Further elaboration of the theory will, it is 
believed, show similar agreement at higher temperatures. 

At 17°C. the rate of hatching is so low that it can be neglected, and, as 
Figs. 2 and 3 show, the population curve follows very closely the theoretical 
curve of Fig. 1. It should be noted that the theoretical curve of Fig. 3 was 
not drawn from values computed from Equation (29), as, owing to lack of 
knowledge of several biotic constants, it is impossible to obtain such numerical 
solutions. It was however drawn on a basis of the information obtained as to 
its shape, from Equation (29). The claim to agreement is thus on the form of 
the curve only. 


The Population Growth Subsequent to T=?;(1) 


Subsequent to T=4,(1), z.e., after the first hatching of eggs commences, the 
problem assumes much greater complexity, due not only to the presence of 
larvae which act as predators, but also owing to the complicated functions 
descriptive of the rates of transformation of each life-history stage into the 
succeeding one. 


The general form of the equations 


By an extension of the reasoning used to determine Equation 18, it will be 
seen that the rate of change of the egg population is equal to the basic rate of 
laying of eggs, namely, ReNy,, minus the rate at which eggs as individuals are 
being eaten by the adults (C’s92+C’s;2), minus the rate at which eggs are 
being eaten by the various larvae of different instars (C’;2,C’s2,:... C's), 
minus the rate at which eggs are arriving at hatching age (Z;)._ It should be 
noted that not all eggs which arrive at hatching age actually hatch. A certain 
percentage shrivel up and die at this time. This matter will be referred to 
later. 

Thus, in a purely general way, where all six larval instars are present 

i=8 
Se « Raliy~ 2 Cpa Oag- Caged (34) 


i=m3 





Certain assumptions are necessary in order that the differential equations 
descriptive of the growths of larval populations may be written. Some 
changes do occur in the biotic characteristics of a larva of any given instar 
during the time that it spends in that instar but, as these changes are relatively 
small, it will be assumed that the various parameters descriptive of each instar 
are constant throughout the life of that instar, and that the change from instar 
to instar is abrupt and discontinuous. 

It has been stated that not all eggs arriving at hatching age actually hatch. 
The same thing applies to all subsequent transformations. Hence if U3, U,.... 
Un be the percentages of deaths at egg-first instar, first instar-second instar, 
etc., transformations, then, when a form (i—1) is arriving at the age of trans- 
formation to a form i at a rate Z;, the form 7 is coming into existence at a rate 
U;Z;j. 
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Then if N, be the number of larvae of the form having subscript j at any 
time 7, and C’,;, be the rate at which they are being consumed by other 
larvae of forms having the general subscript 7, etc., then, in a purely general 
way, 





aN; = 
= = U;2Z; ~2 C'5.5- C10. 5 -— C1 5 — 2541 - (35) 


It will be seen later that some or all of the C’’s may be zero, depending upon 
the circumstances at a particular time. 

It is apparent that an equation similar to that descriptive of the growth of 
the egg population can be written in the case of the growth of the pupal 
population. 


Evaluation of the C’’s 


It is assumed, as again is practically true, except in the case of some of the 
smaller larvae, that no larva of a given instar can consume one of the succeed- 
ing instar, and that adults are immune from attack by any form, even by other 


TABLE VIII 





es Ds “gs Mee Ses ae Me ee Bee 
Quis Qi1.2 Qus Que Qus Qi1.6 Qi1.7 Qiis Qi1.9 
X's X'so2 X's Xie Xs X'sos X'or X's Xo 1 1 
Qi0.1 Qr0.2 Qio.s Qro.4 Qio.s Qr0.6 Qio7 Qi0.8 Qr0.9 
1 1 i 1 1 1 1 1 1 1 1 








X's X's2 X's.3 X's X'ss X's6 X's7 X's X’sn «611 
Qs: Qs.2 Qas Qs4 Qs.s Qs.6 Qs.7 Qs.s Qs.9 














Ts Te fs Eu 2s Su Bu 1 Ry 1 1 
Qr1 Qr.2 Qr3 Qr4 Qrs Qr6 Qr.1 Qr9 
tk Ca te fe ee +. + mh + 
Q6.1 Qe.2 Q6.3 Qs.4 Qe.s Q6.6 06.9 
Ee ty fe Eu ke 44.04 ga 
Qs.1 Qs.2 Qs.3 Qs4 Qs. Qs.9 
a Oa a 1 1 ee ok 
Qa Qs2 Qas Que Qs 
eh” Of 1 1 1 a a 
Qs Qs.2 Qs. Qs.9 
i i i 1 1 1 1 1 1 1 1 
1 1 i 1 i 1 1 i 1 1 1 


* It seems doubtful if these are of much importance since Ps 5, Pg 2nd P34 are small and prob- 
ably sero. 
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adults. (This latter assumption is strictly true, except for a very brief period 
during which a new adult is emerging from the pupal skin.) These assumptions 
may be introduced into the formulas by writing P;;=0, where j>i, (except 
j=9), and P; 40, Pi1;=0 for all values of 7. It is also evident that, of necessity, 
as flour, eggs and pupae cannot eat, all P’s of the forms P,,, P,,; and Py; 
are zero. 

The C’s may now be evaluated by an extension of the reasoning used in 


X's. : 
Equations 19 to 24. As such symbols as “— occur frequently, they will 
8.3* ij 





be written as! - 
4.§ 
Certain of these are of no importance, namely, those for which P;;=0. 
Such a case gives, of course, the indeterminate form -<., but it is easy to see that 
this reduces to the value 1. There result, then, the ratios shown in Table 
VIII, all equal. Those for which P; ;=0 have been written as 1 immediately, 
and thus disposed of. 


As, (See Equation 20) 
jwii* 


EN,= X X'5, 
j= 


and as any one of the X's may be obtained in terms of any other in the same 
row by means of the relationship 


X54Qi. . 








XG = (36) 
Qin 
ite Xn 
on substitution, EN;=—= D Q;, (37) 
ik jot 
ENQis | 
whence X’, = Jail (38) 
” = Qi.5 
jut 
which, on dividing by W,A,, becomes 
ENQis 
Cia = — (39) 


WA, X Qs 
j=l 


It is now necessary to consider the forms of the Q’s. Some are equal to 
zero directly, and are thus of no further interest. These are Q’s of the forms 
Qi.j, Q2.4 Qos, Or, where i<j (except j=9), and Q;10, Qj: for all values of 4. 
It must be remembered, of course, that any Q; ; is zero if N; is zero. 

As a preliminary step, the Q’s may be divided into four classes, (a), (b) and 
(c) as on page 642, with an additional class (a,) for encounters with flour. Thus 
remembering that Q5=P;; Y's 


*The summation may be made from 1 to 11 since a and Xi107=9 ‘i 
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Under case (a), 


P;,A,;W,Kguser5 1 N3G 
Go Seeet eee (40) 


PsA WiKimeiaNG mri NG ‘ (41) 


1 = G { 


Puy 14 Wi Kyu trina NiG . 


Q141=— G (42) 





Under case (a) there are formulas descriptive of amounts of assimilable 
material obtained by larvae or adults from eggs or pupae, thus 


2 
P32A2W2N3N2u3573.2 


= (43) 
Q3.2 S 
Qs = P42A aM aMelorartas te (44) 


2 
Pi9.242W2NioN2109*"10.2 : 


Q10.2= (45) 
10.2 G 
PocAgWoNoNousere 
0s. sa 6.94 gWoNoNoustte.9 : (46) 
G 
P11 9AoWoNi1 Noms =? 
One™ - 11.949 oNiiNomii F119 (47) 


G 


Under case (b) there are formulas descriptive of amounts of assimilable 
material obtained from larvae of a given age by larvae of the same age, thus: 


4P; 54 3W3N3(N3—1)us*3.3 : 


a 48 

Q3.3 3G (48) 
4Py AgW iN No—1) ugar 

— Pe ee ; (49) 


7 4P; ,A;WNA(N,—\) meri 
3G 


- 4PsAgWeNa(No—1)usr758 ' 
3G 
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Under case (c) there are formulas descriptive of amounts of assimilable 
materials obtained from larvae of various ages, by larvae not of the same ages, 
and by adults, either mature or immature, thus: 

Py 3A3W3NaN3F (ug us) *743 


Qs = ee ee , (52) 


nm mess us) "753.3 (53) 


P WoNiurNoF(ury, us) re 
Case 11.64 6WoNiiNoF(urs, M6) 11.6 (54) 


G 
The C’’s may now be evaluated in terms of the Q’s. The formulations may 
be simplified as a result of an examination of F(u;, 4;). It has been convenient 
hitherto to use the separate velocity corrections; it will now be shown that they 
may all be written in the form F(y,, 14;). 


[ —dae—mb? +m? (5S) 





F( uj, u;) = 
6ujuy 

Suppose that 4; >; 
Then if 4; or 4; is zero, while the other is not, clearly it must be u; =0. 


ae 
Whence Pi. a) = Soe (56) 
Fujuy 


Differentiating both numerator and denominator separately with respect to py 
we obtain, 


3urt+3uy ; 
3m; 


Whence uj Fis 1g) = mye 


Similarly it can be shown that where y; > 4;, 
ki 
so Fy my) = aye 
It is also evident from direct substitution in (55) that, 
hi 4 4 
wu, Flnin ay) = 3 Mi (Or > a) 
Whence, where i1=3, 4, 5,....8, 10 and 11, 
Py ,WA,K,F(uj0) 97, .NGG E, 
Cus S — 


WA.) PiaWrArKsPluO) er NG 7S) Py WiAgP(wintsy) 001 3NGNy 
as: Tneremenney meee te ecare ad 
j=2 
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Pi 1K sF( uj, 0) 4.1 ENG 





al 2 - 2 (58) 
Pi Wi A, K,F(uj, 075 5G+ LD Py sWjAj Flu, my) iN; 
j=2 
and where j=2, 3, 4, ....9 and j¥i 
- P;5F (ui my)02 sNNjE, 
Ci.3= 


kad (59) 
Pi) WyA,K,F (uj, 0) 134G+ = PinWeAn Fis unten 


and where j =1=3, 4, 5, 6, 7, 8, (C’;,; =O when i=1, 2, 9, or j7=10 or 11) 


Py F (ui, mire; N; (N;—-1)E; 
kad (60) 


Pi WyAyF(uy Ores G+ SS PypWeAn Flug, unten Np 
k=2 


C7147 


from which it is readily seen that C’; ; is zero if either N; or N; is zero. 


Determination of the forms of the Z’s 

As stated above, Z; is the rate at which a form with subscript (¢—1) is 
transforming to a form (4) at a time 7. 

Since the Z’s involve certain functions the forms of which have not yet been 
determined, only a cursory explanation of the forms of the Z’s will be given 
here, a more detailed explanation being reserved for later publication. 

By an extension of notation previously given*, let fe, ts...... ti, be the fixed 
times at which transformation to the forms having subscripts 2, 3, 4,...... il 
commences, t.é., the first pupa comes into existence att. Alsolet I, I3,...... 
I be the durations of individuals in the various life-history stages having 
subscripts from 2 to 10.f 

It will be clear that the history of an adult, for example, reaching maturity 


10 

at a time T may be summarized as follows: laid as an egg at T— 4 J';; hatched 
: 2 
10 

to a first instar larva at T— SI; transformed to a second, third, fourth, fifth 
3 


10 10 10 10 
and sixth instar larva at the times, T— D7, T- DI, T-DI,T- Di 
4 5 6 7 


10 
T— DTI, respectively; pupated at T—I'—Iy; emerged as an immature adult 
8 


at T—I'y; and became mature at J. Thus with each individual there are 
associated past times at which its various transformations took place. 
Consider two groups of eggs, a and §, one of which, a, hatches promptly 
at 4, while the other hatches later. In the case of the first group, a, the de- 
predations due to eating have been caused only by the Ny adults originally 
present at T =,=0, whereas in the case of the group 8, in the early part of its 


*As only first generation eggs will be discussed, we may write te, ty, t4, rather than t,(1) t(2) etc. 
tI'n ts indefinite, being the life of the mature adult, as such. 
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life as eggs it was attacked only by the N,, mature adults, while in the later 
part of its life it was attacked by both the N,, mature adults and the N; first 
instar larvae existing subsequent to the time 4;. It is thus apparent that any 
function embodying the rate at 

TABLE IX 


which this group, 8, has been eaten VALUES OF /2 TO fy) AND I; TO I‘) AT 27° C. 
at various times, must be discon- ; a 























tinuous at fs. If moreover, T be ¥ Value | t | Value Derivation of ¢ 
greater than t,, discontinuities  ——— 

must ‘also exist at such of the Tr |6.04| & | 0.00] h=%=0 

times Se are ty, as fall within Ts 2.43 ts 6.04 | Atl: 

the time span over which the group r He : ae ores ~~ 
has been in existence as a group To | 3.27] & | 15.13 | a+0a+ rs 
of eggs. Similar discontinuitieswill 7 |é@;/ 2 | 2s letnt 7 
occur in other functions embody- Ty | 8.64} t | 28.46 | a+Iit+..... I, 
ing the rates of eating of other Pr or 2 ae et Sas rn. 
forms, and since the Z’s are such | 


functions, the above discontinuities 
will occur in the Z’s as will be shown below. For reference the values of & 
to t, and of I’; to I are given in Table IX, computed from Tables I and V. 
It has been shown (Equation 31) that 
lim 4N; 


T—h=0 “aT =ReNu, 


that is, eggs are coming into existence at this rate when T=4,=0. 

Consider the rate at which eggs are arriving at hatching age when T=%,. 
Clearly this must be less than ReN;, since many eggs have been eaten. 

Let this rate be O(¢s). 

Then the total reduction from the rate at which eggs were coming into 
existence at T=#=0 to the rate at which they are arriving at hatching age 
at ft; is ReNy —O(ts). 

But the rate of eating of a small group of eggs laid between T =t and T =i+dt 
is, at any time ¢, 

HOC 2 4 


N; 


whence 
ty 


a(t)C’ 
Ot) =ReNy — [ *Aente (61) 
N2 


When T=4;, first instar larvae (subscript 3) are present, whence during a 
sufficiently small neighborhood where T>t, 
T 


by 
6(t)C’ - - 
ur) =ReNs— | OC a [ OCustOr) mn 


N, N; 
T-Ts b 


From Table IX it will be seen that 4:=6.04 and that 4=8.47, whence, 
as I;=6.04 and J';=2.43, it is possible for T to exceed & without T—I; 
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exceeding ¢;. In such a case, fourth instar larvae are present during the perioa 


from & to T and we have 
ty he 


f a(tC 2 a(t) (C’ ~ 
6(T) = ReN,,— | Soll de x | 0 (C2 C's2) 
J N2 N2 
1-Ts ty 
y 
~: | OO)(CrratCu2tC'a2) 4 
N2 


dt 


(63) 
t 
When 7T=12.08, T—I°:=¢% and the first integral on the right hand of 
Equation 63 becomes zero, since the upper and lower limits are equal. 
When T becomes greater than 12.08, 7—J°; becomes greater than #; so that 
the integral becomes negative. Since this has no connection with any real 
case in an actual population, the fact that only positive integrals are to be 


+ 
considered will be denoted by a plus sign, thus: | , 
As for a sufficiently small neighborhood, T> 12.08 


te 7 
+ " ’ + ’ - ’ | 
(7) =ReNy, - / HO (CustCos) 4 / ONE ngt Carters) a, ae | 
T-Th , ta , 


It is seen that as T increases, there is a periodic addition of integrals on the 
right hand, and a concomitant extinction of integrals on the left, accompanied 
by a migration of the lower limit T—I,._ This migration will be denoted by an 


+ 
arrow, thus: 


T—Ts 
nei 


Whence, for a value of T greater than say ¢; where ¢t,; is any one of the ?’s 


ts 
+ v +/ ’ 
0(T) =ReNy,— / OC 2 9 | a(t) (Cn2tC's2) 4, 


N; N: 
T-—Ts ts 
bs jus T j= 
*f o(t) (cr +22) +f o(¢) (Cus +3c2) 
= acide Decals. 6. im Ri Hi sist Ranieasale asia dt, (65) 
be Nz 4 Nz 
which may be written in contracted form as 
a ‘ + ett j=uk 
0(T) =ReNy, - / 6(t) Cu. yy > 7 (cus#3cia) 
Ny k= 3 dt 
T-Ts bb Nz 
—> 
+ T j= 
08) (Cnet 
F ( ( 11.2 Z 2) " (66) 


‘ Ne 
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It will be apparent that similar functions may be written for the rates of 
transformation from first instar to second, etc., but as little is known as yet 
as to the form of the function 6(t) these formulations will not be included here. 


The egg and larval populations for a sufficiently small neighborhood, T >t; 

When an egg hatches, the larva is at first comparatively inactive, and as the 
only material lost during the transformation from egg to-first instar larva is the 
minute amount which makes up the shell, during a sufficiently small neighbor- 
hood,’ T >#s, first instar larvae may be thought of as eggs. In this case it will 
at once be apparent that the sum of the two populations, eggs and first instar 
larvae, will be equal to what the egg population alone would have been, in 
the absence of hatching, with the exception of those eggs which entirely pass 
out of existence owing to natural mortality at hatching. That is to say, if 
Nz be the value which the egg population would have reached at a time T, in 
the absence of hatching,.and if Ne and N; are the actual numbers of eggs and 
first instar larvae, at any such time, 


N,-—(1—U3)Z3=N2+N3. (67) 


If now, as is the case, (1— U3) be small, for a short time after hatching com- 
mences, two situations may arise: 
(a) The egg population has reached a point of quasi-equilibrium, 7.e., is 
oscillating between the greatest integral number less than £ and the least 
integer greater than ¢ (see Equation 30). In this case as N’;>0 when T>4;,N; 
will become negative as soon as T >4s. 
(b) Such a condition of quasi-equilibrium has not been reached by the egg 
population when T=%. In this case, the egg population will continue to 
increase until such a value is reached that 
ink 
ReNy=Z3+BCiat Cur (68) 
where k is the subscript for the oldest larva present at the time that the value 
of Nz in Equation 68 is reached. After this time, the egg population decreases. 
Inactivity of first instar larvae is not necessary to the truth of the above 
expression however. 
It seems hardly necessary to point out that immediately subsequent to 
T =ts, the expression 
i=k 
Usds> BC, +C'113 


must be satisfied if the larval population is ever to come into existence at all. 
The question might now be asked as to what would occur should trans- 

formation to second instar larvae, (subscript 4) not take place. So far the 

writer has not been able to determine this point mathematically, but it would 

seem that there are two possibilities: 

(c) The egg and larval populations come into stable equilibrium, in which 

the rate of egg transformation to first instar larvae U;(Z;3) at all times balances 
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the larval loss due to natural death, and due to the eating of larvae by larvae 
and adults, and in which Z;, and the eating of eggs, and the death of eggs are 
together equal to ReN,,. That is, N; and N, are zero. 

(d) The egg and first instar larval populations enter into an infinite series 
of ‘‘cyclical’’ changes, where by “‘cyclical’’ it is not necessarily meant that 
each cycle is identical with the preceding one, but is nevertheless of the same 
general form. 

It is felt, lacking rigorous mathematical proof, that (c) can be ruled out for 
the following reason. At any time TJ the rate of production of first instar 
larvae is a function of all the past history of the egg and larval populations 
during a period of time from T—T, to T, that is to say, 


a 
+ 
, "gs 
vsts=0s| Rev | ay (Le2t us Ja], etc. (69) 
2 


T-T: 
Suppose that N2 and N; have just arrived at such a pair of values that 
U3Z3=C'33+C'13, (70) 


4.e., the equilibrium mentioned under case (a) is attained. Let T receive a 
finite increment AT, then, 


T+AT 
"dg , 
Us2y=Us| Re / 6(t) (£s24 Cus Ja | (71) 


2 
T+AT-T? 


Then if UZ; still equals C’; ;+C’s,;, Ne and Ns must have had the values 
mentioned above throughout the interval T—I, to 7, whence the above 
equilibrium cannot be maintained for a finite time unless it has already been 
maintained for a period I>. Since Nz and N; are not at these equilibrium 
values when T=%;, they can never, on attaining them, remain at them for a 
finite time. It follows that where Z,=0, i.e., where no transformations to 
second instar larvae occur, the first instar larvae increase in numbers until 


U3Z3=C'33+C'13, 


and then, due to eating, decrease again. A second such maximum cannot 
occur again until a period of at least J, has passed, for not until then can the 
increase in egg population (resulting from a decrease in C’;,) have any major 
effect upon the larval population. 

It is true, of course, that there will be a retardation in the rate of decrease of 
N; due to the increase in Z; with decrease in C’;, but this cannot in general 
make N;>0. 

Generally speaking, where Z,>0, the above-mentioned maximum value of 
the first instar larvae is not attained, owing to the reduction caused by Z,. 

It is believed that the other larval instars increase in numbers, tending 
towards analogous maxima, and in general, decrease in numbers before reaching 
them, for analogous reasons. 
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Distribution by ages in the total larval population 

During actual counts of larval populations which had been running for some 
time, it was noticed that there was an increase in the proportion of older larvae, 
even after the time when larvae of all stages were present. That is to say, 
with older populations the distribution by ages tends to become increasingly 
skew in the direction of the greatest age. The explanation of this phenomenon 
from the theory thus far developed is comparatively simple. 

Consider a group of first instar larvae hatching from eggs during a period 
ts< T<%+At, where At is some small increment of time. By reason of the 
postulate that C’, ,=0, i<j except j =9, all through their lives this small group 
is subjected to an “eating force” of only C’;;+C’; ;, where j is the subscript 
indicative of the instar in which the group exists at any time. On the other 
hand, a similar group hatching out over a period of equal duration, 
but situated later in time, when there are perhaps larvae of the second 
and third instars present, will be subject to much greater eating, namely, 
C’ 425 tC GengtC 55+C' 15. Not only this, but the group will be 
smaller to start with, since jt is easily seen that in populations of the type 
oT 
Thus, theoretically at least, the first larva hatching promptly at T= will 
always pupate, provided it is not eaten by an adult, and does not die from some 
natural cause. Obviously then, with the passage of time, all but the older 
larvae tend to disappear, and the age distribution is skewed in the direction 
of the greatest age. 





described in this paper <0 where 2 is the subscript of any larval instar. 


Pupal populations 

In general it may be said that the pupal populations tend, as do the larval 
populations, to reach a maximum value, at which the increase due to trans- 
formation from sixth instar larvae (subscript 8) just balances the loss by natural 
death and by eating, but do not attain this value by reason of transformation 
to the adult form. 


Adult populations, T >to 

The equation descriptive of the growth of the adult population is in general 
far simpler than those descriptive of the growths of the egg, larval and pupal 
populations owing to the fact that C’; 10, C;4, =0. 

Therefore, the only negative term will be one to account for the slow and 
gradual death of adults by reason of sheer accident and infirmity. This term 
may be written in the form K (7—#) where K is a constant. 

What then limits the size of the adult population? The equation may be 
written as 

Nyy =Zyy—K(T—t1), (72) 


whence it is at once seen that the limiting factor, since K is very small, is a 
diminution in Z,,, the rate of transformation from immature to mature adults. 
Thus, if K be considered as practically negligible, Nu will increase until 
Z1=0 which practically, as C’; =0, means Z;)=0. 
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An interesting point arises here. Suppose that Y adults, mature and 
immature, are sufficient to reduce and hold the pupal population practically at 
zero, so that N’ is zero. By reason of the fact that some time is necessary to 
eliminate the pupal population in existence at 7 =f the adult population 
generally rises above the value Y which would obtain under any particular 
set of environmental conditions (temperature, flour mass, etc.) As this large 
adult population cannot reduce the pupal population to less than zero, it feeds 
upon flour, eggs, and such larvae as emerge from time to time, and in the mean- 
time dwindles slowly away (see Figs. 4 and 5) by reason of the second term 
in Equation 72, to approach the value Y from above. 

Owing to unavoidable failure in the mechanical equipment used to maintain 
constant temperatures, no populations were grown for a sufficient length of 
time to determine what would happen beyond this point, but in the case of 
some populations in which the number of adults was accidentally greatly 
reduced, a new cycle analogous with those mentioned on page 660 commenced, 
the adult population rising again to a value in excess of Y. It is believed that 


: ; ‘ Nu-Y os 
a periodic fluctuation with a period of t,+ = would be set up, where Nui 


is the maximum number of adults reached under any particular set of condi- 
tions, and ¢, is the time from & to the time of reaching the maximum number 
of adults. 

It might also be asked how hatching of eggs can fail to occur, as it at times 
does, when subsequent to T=. there are perhaps 1000 eggs present. The 
answer to this lies in the fact that Z; is dependent not only upon the number of 
eggs, but also on the rates of consumption of the same throughout their life as 
eggs. Ina case such as the above, the rates of consumption are such that the 
probability that an egg will be found and eaten in a time less than 7; is 1. 
Hence no eggs survive long enough to hatch, even though they are very numer- 
ous. Such a condition is of course unstable, in that a diminution of larvae or 
other predators will cause an immediate increase in the numbers of first 
instar larvae. 


The Protective Influence of One Form on Another, in the Matter of Eating 


Consider the addition of a new form to the system at any time. It will be 
apparent that if this new form is edible, it must bear some of the depredations 
of the predators, and for that reason must diminish to some extent the rate at 
which all other edible forms are consumed. The mechanism of this, mathe- 
matically, is quite simple, as it operates through the fact that the addition of 
a new edible form adds another Q to the denominator of all the C’’s already in 
operation, and since the added quantity is real and positive, the C’’s are 
thereby diminished. If it should occur that the added form is predatory as 
well as edible, the protection offered may be more than offset by the additional 
predatory influence. This will of course depend on the relative value of the 
new C’ introduced and the diminution in those already present. 

This matter will be referred to more extensively in Part 3. 


cea 
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Part 3 
Discussion of Experimental Data 


In Tables X, XI, XII and XIII are tabulated the actual population counts 
from duplicate experiments at four temperatures, namely, 17°, 22°, 27°, and 
32°C. all at 75% relative humidity. (This was somewhat exceeded in the 
cases of the two populations at the lower temperatures.owing to unavoidable 
condensation on the cooling coils in the controlled cabinets.) The duplicate 
experiments are denoted as (A) and (B). The mean values from (A) and (B) 
are also shown. 

These means are plotted as Figs. 2, 4, 6 and 8, while Figs. 3, 5, 7 and 9 show 
curves drawn on a basis of the actual data, but with modifications made in an 
endeavor to replace information lost by reason of the relative infrequency of 
the counts. These changes have been made on the basis of the information 
gained from the theory developed in Part 2. Care has been taken to make 


TABLE X* 


POPULATION SERIES IN WHOLE WHEAT FLOUR, 32 GM., AT 17°C. AND 75% 
RELATIVE HUMIDITY 
































| | 
Time, | Eggs Eggs | Eggs | Larv. | Larv. Larv. Ad. Ad. | Ad. Tot. Tot. Tot. 
days | (A) | (B) | Mean (A) | (B) | Mean | (A) | (B) | Mean | (A) | (B) | Mean 
oO 16 16 16 16 16 16 16 
10 41 24 | 37 16 16 16 57 | 40 48 
20 61 47 54 | | 16 16 16 77 | 63 70 
30 102 93 97 | | 16 16 16 118 109 113 
40 134 112 123 16 16 16 | 150 128 139 
50 162 141 152 16 16 16 178 157 167 
oo 186 | 160 173 2 1 16 16 16 204 176 190 
70 178 153 165 1 a 16 16 16 195 168 182 
80 219 143 181 3 2 2.5 16 | 16 16 248 161 204 
90 245 130, | 187 | 5 0 2.5 16 15 15.5 266 145 205 
100 253 149 201 5 3 4 16 15 15.5 274 167 220 
110 249 168 208 4 0 2 16 15 15.5 269 183 226 
120 214 169 191 13 1 7 16 15 15.5 243 185 191 
130 222 168 195 13 1 7 16 15 15.5 251 184 217 
140 309 189 249 14 0 7 15 is | 45 338 204 271 
150 275 203 239 12 0 6 15 1s 15 302 208 255 
160 190 186 188 8 G 4 15 15 15 213 201 207 
170 212 213 212 8 5 6 15 15 15 235 243 239 
180 203 220 211 6 5 5 15 15 15 224 240 232 
190 243 265 254 8 3 5 3 | #138 15 266 283 274 
200 258 230 244 15 9 12 1S | 15 15 288 254 271 
210 219 245 232 16 19 17 15 1s 15 250 269 259 
220 210 224 217 31 39 35 15 15 15 256 278 267 
230 261 233 247 32 32 32 15 15 15 308 280 | 294 
240 221 201 211 8 57 32 | 15 15 15 244 273 258 
250 229 190 209 20 82 1 +s 15 15 304 287 295 
260t | ... e. i ei a Ie, Eis a at wa Efe a 
270 234 174 204 42 138 90 15 15 1S 291 327 309 
280 219 149 184 40 95 67 | 15 15 15 274 259 266 








* The writer is indebted to Dr. R. N. Chapman for the data in this table. 


+t Count missed on this day. 
NotE:—No pupation occurred at 17°C. 
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only such modifications as could be found homologously at all temperatures 
with the exception of 17°C., in which case, owing to the great value of I; 
only the earliest stages of the population growth are manifest. 

As soon as T becomes greater than #,=0, eggs commence to accumulate 
according to Equation 29, and as has been shown, tend towards the limiting 
value £, at which point they are found and eaten as rapidly as they are laid. 
At 17°C., owing to the high value of J2, this condition is practically attained, 
but as at times, a few eggs do hatch, in the region of 60 to 70 days, Fig. 3, 


N’2 = REN, —Z3—C'3.2—C'1,.2. <0. (73) 


At 22°, 27°, and 32°C. this equilibrium value £ is not reached owing to the 
reduced value of I2,and at p, the egg population reaches a maximum as described 
in Equation 68, and then decreases owing to hatching to form first instar 
larvae. 

Thus, the population at 17°C. is an example of Case (a) page 659 while the 
other three populations illustrate Case (b), page 659. 


TABLE XI* 


POPULATION SERIES IN WHOLE WHEAT FLOUR, 32 GM., AT 22°C. AND 75% 
RELATIVE HUMIDITY 
























































Nl l l 

Time, | Eggs | Eggs tiinkiae’ Larv. bi Pup. Mage Pup. | Ad. | Ad. | Ad. | Tot. | Tot.| Tot. 
days | (A) | (B) | Mean| (A) | (B) |Mean! (A) | (B) | Mean| (A) (B) |Mean| (A) |e. (B) |Mean 

Pape ean ee nee gern gery pe ape ie omni end 
0 | | | 16| 16] 16| 16| 16] 16 
10 | 196 | 222} 209; 11| 12! 41 16| 16| 16] 223 | 250} 236 
20} 338 | 314| 326| 79| 661} 70 | 16| 16 | 16| 433 | 391 | 414 
30 | 399] 383! 391 | 210| 131 | 170 16| 16| 16] 625 | 530| 577 
40 | 371 | 398| 384) 323 | 219 | 271 0} oO 0| 16] 16| 16] 710 | 633| 671 
so | 287 | 281 | 284] 315 | 204 | 259 i .8 2] 16] 16] 16] 619] sos! S61 
60 | 225 | 186 | 20S | 365 | 326 | 345 7) 9 8| 16} 19] 17] 613 | 440] 526 
70| 200| 149 | 174] 372} 248/310; 22} 22] 22] 19] 26] 22] 613 | 445] 529 
80 | 132 118 | 125 | 297] 153] 225] 97] 45| 71 | 27| 35] 31] 553 | 356] 454 
90 | 79) 124) 101/175} 176) 175| 93) 81 | 87] 76) $9] 62} 423) 440] 431 
100 | 106 | 159} 132| 92| 129/110) s7| 74| 65 lisa] 83] 98! 369 | 445 | 407 
110 | 161 | 243 | 292 | 47| 82/| 66 | 10 | 58/34 | 136} 124] 130| 354} s07| 430 
120 | 157| 351; 254| 45| 83! 64! 2 | 37 | 19/130} 142 | 136 | 334 | 613 | 473 
130 | 128 | 347 | 237| 69/| 176|122| 2] 9 S| 120] 148 | 134 | 319 | 680 | 499 
140 | 151 | 325| 238 | 72| 131 | 101 o| 2 1 | 109 | 140 | 124 | 325 | 598 | 461 
1s0| 148 | 235 | 191 | 77] 185 | 116 o| oO © | 101 | 137] 119 | 326 | 527| 426 
160 | 87 | 158 | 122] 67| 162) 114| 0] 0 0| 92 | 128 | 110| 246 | 448 | 347 
170 | 250 | 382 | 316] 37 167 | 102 0} oO 0} 87| 120| 103 | 374 | 669 | sai 
180} 621 | 610} 615 | s2{ 122} 87| 0 | 0}; o| 86| 119| 102 | 759 | 851 | 80s 
190 | 487 | 336 | 411 | 207 | 201 | 206; 1| 1 1} 80 | 121 | 100! 774| 659| 716 
200 | 377 | 224 | 300 1248 | 182/215; O| 4 1} 76| 118| 97] 701 | 525 | 603 
210 | 224 | 134 | 179| 346/ 201/273 | 5 2| 71 | 114] 92] 641 | 454] 547 
220 | 144] 65 | 104| 328; 184 | 256 o| 9 4| 69| 106| 87] 541 | 364| 447 
230| $3) 49] 66 | 307| 148 | 227 6 | 2 4| 63} 88] 75| 409] 287] 348 
240} SO} 37| 43{279] 96] 187/ 13] O 6| 62} 80] 71] 404/213! 308 
250} s4| 34] 44| 263| 59! 161 7| 1 4| 60] 64] 62] 384| 158] 271 
260 | 268|158| 218/158 | 23] 89 0} o 0} s6{ s5| sS| 479] 236) 357 
270 | 284 | 554| 419 | 230] 44! 137 o| o o| ss| 48| 51| 669 | 646| 657 
o| o 0| 55| 46] so| 876 | 828| 8s2 


280 | S21 | 580 | 550 | 300 mie 
| 


* The writer ts indebted to Dr. R. N. Chapman for the data in this table. 
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The writer confesses inability to explain the sudden fall in larval population 
at p: in the case of the population grown at 22°C. Lacking knowledge of the 
value of the J’’s for 22°C., it is not possible to determine whether this is an 
inherent peculiarity of the population growth or whether it is due to faulty 
technique. It is possible that it is due to relatively great activity on the part 
of sixth instar larvae at this temperature, but this is hardly in accord with 
observations made of such larvae. Moreover, such activity would manifest 
itself in an even sharper decline in egg population. Again, had such a decrease 
in the larval population actually occurred, it would have been reflected in a 
similar drop in pupal population at a time J’; later. As there is no evidence 
of this, beyond the fact that the pupal curve seems somewhat bluntly truncated 
at 90 days, the time of the expected discrepancy, it is assumed that it is due to 
an error in count of an even hundred in ‘‘B”’ at 50 days, the finely waved line 
indicating the appropriate correction. 

T now increases to approach and exceed fs, and as pupae make their appear- 


ance, 
i=8 


N'9=U9Z9—LC'i9—C'11,9> 0. (74) 
in3 
TABLE XII* 


POPULATION SERIES IN WHOLE WHEAT FLOUR, 32 GM., AT 27°C. AND 75% 
RELATIVE HUMIDITY 



































Time, | Eggs |Eggs| Eggs |Larv.| Larv. |Larv.| Pup. | Pup.| Pup. | Ad. | Ad. | Ad. | Tot. | Tot.} Tot. 

days (A) | (B) | Mean} (A) | (B) |Mean! (A) | (B) | Mean; (A) | (B) {Mean} (A) | (B) | Mean 
| | 

| Po bgt ae Sa 

0 16 16 | 16 16 | 16 16 

10 | 363 | 383 | 373] 18 S1| 34 15}. 17] 16} 396) 451 | 423 

20 | 350 | 299 324| 319 316 | 317 Sb Sib ee 17} 15 | 669 | 633 | 651 

25 | 213 | 227 |. 220 | 377 | 388 | 382 2 6 4} 13 16| 14! 605 | 637/| 621 

30 | 149 | 143} 146 | 365 | 403 | 384 7| 14 10} 13 17| 15} 534] 577] sss 

35 | 259 | 261 | 260 | 309 | 308 | 308 56 | 73 64) 19/ 24 | 21 | 643 | 666} 654 

40 | 202 | 225 213 | 261 | 288 | 274| 171 171 | 171] 26 46 | 36 | 660 | 730} 695 

46 | 331 | 229] 280/ 157 | 183/170 | 138 | 133 | 135 | 161 | 173 | 167} 779} 718| 748 

50 | 341 | 218 | 279 | 128 | 159) 143) 69 | 90 | 79 | 229) 236 | 232 | 767 | 703 | 735 

55 | 334 | 332 | 333 | 144] 142 | 143 26 | 27 26 | 272 | 292 | 282 | 776! 793 | 784 

60} 328 | 331 | 329 | 104| 137] 120) 118)| 2) 19 | 285 | 304 | 294 | 735 | 793 | 789 

70 | 428 | 327; 377) 41 87 | 64 24 | 39 31 | 299 | 333 | 316 | 792 786 | 789 

80 | S22 | 492} 507} 12| 30] 21 17,| 27 22 | 308 | 356 | 332 | 859 905 | 882 

90 | 706 | 869} 782, 10 3 6 3| 27 15 | 303 | 365 | 334 | 1022 |1264 | 1143 

100 | 824 |1040 | 932] 49 36 | 42 4 1} 2 | 317 | 371 | 344} 1194 [1448 | 1321 

110 | 812 | 750} 781 | 73) 109 | 91 a}. 4 1 | 297 | 367 | 332 | 1183 |1227 | 1305 


120 | 509} 615 562 | 221 157 | 189 15 1s 15 | 316 | 374 | 345 | 961 |1097 | 1029 
130 | 388 | 460 | 424 107 100 | 103 27 40 33 | 321 384 | 352 843 7 984 | 913 
140 | 491 | 717 604 | 92 64 78 16 | 35 25 | 325 390 | 357 | 924 |1206 | 1065 
150 | 323 | 681 502 7s 47 61 18 14 16 | 323 405 | 364 739 |1147 943 
160 | 216 | 274] 245 73 98; 85, 14 | 8 11 | 332 416 | 374 | 635 | 796 715 
150 | 372 | 446| 409 | so 74 | 66 2 12 7! 326 | 412 | 369 759 | 944 851 
180 138 | 110 124; 30 57 43 0} 13 6 | 319 | 401 | 360; 487 | S81 534 


1904 























* The writer is indebted to Dr. R. N. Chapman for the data in this table. 
+ Population accidentally killed at 190 days by escape of ammonia from refrigeration equipment. 
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This pupal population tends to increase to such a point that, 
i=8 
UsZ9=LC'i9—C'r1, (75) 
i=3 
but by reason of transformation to immature adults, arrives only at the 
maximum at p; where 


ims 
Ustg=LC'sa—C'11.9—C'r09—Zi0. (76) 
The larval population continues to increase subsequent to T =, for a short 
period, reaching a maximum at 3; where 


dull im8 
N’,=U3Z3-2,-L L C,,=0* (77) 
imj jus 
Subsequent to T = and also to the time of 3, the egg population continues 
to fall, to reach a minimum where, 


ims 
N’2=ReNyy—Z3— D Ciia—-C'una- (78) 


i=3 
TABLE XIIlIt 


POPULATION SERIES IN WHOLE WHEAT FLOUR, 32 GM., AT 32°C. AND 75% 
RELATIVE HUMIDITY 


Eggs | Eggs | Eggs |Larv. baile Pup. | Pup.| Pup. | Ad. | Ad. | Ad. | Tot. | Tor 
Mean | (A) (B) |Mean| (A) _— Mean | (A) (B) |Mean/ (A) 


16| 16 
355 | 182 , 1is| 16 
107 | 423 14} 18 
117 | 274 15 14 
149 | 122 89 96 
358 231 | 307 
307 343 | 358 
352 366 
| 357 | 368 
376 | 371 
362 
358 | 371 
361 
353 


13 
925 il 
55 
539 139 
738 106 
906 58 
573 731 37 
280 oo 
614 | 497 $55 | 87 55 
919 |1032 975 | 122 191 
SiS | 190 3S2 | 151 242 
| 506 187 346 | 106 172 
j 
+ The writer is indebted to Dr. R. N. Chapman for the data in this table. 
Note:—Populaiion accidentally destroyed by fatlure of temperature control, at 203 days. 
* Note that Nyw=0 and thatC’s,;=0. 
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and then, by reason of reduction in the numbers of the highly predatory sixth 
instar larvae (subscript 8), and by reason of the protective influence of the 
edible pupae (see page 662), acting as an alternative food for the remaining 
predators, N’, becomes greater than zero, and the egg population again 
increases. 

In the case of the population grown at 22°C. there is an additional refutation 
of the idea that sixth instar larvae grown at 22°C. are highly predatory as the 
diminution of their numbers does not decrease the rate of eating of eggs suffi- 
ciently to cause a great rise in the egg population. Possibly sucha rise did 
occur, and was missed on account of the comparative infrequency of counts. 
It is an unavoidable misfortune that no count was made at 65 days. 
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Fic. 4. Population growth of T. confusum, Duv. at 22°C., mean of (A) and (B) populations. 










T now increases to approach and exceed ti) when immature adults (subscript 
10) begin to come into existence. For a short time subsequent to fy the egg 
population continues to increase, until with increase in Ni a maximum is 


reached at ps where, 
i=8 


N’, =RENyy—Z3—Y C2 —C'0.2-C'11,.2. (79) 
i=3 








Subsequent to this time, for the same reasons, N’’; becomes less than zero, 
and the egg population falls to reach a minimum at -. which will be discussed 
later. 

In the meantime, subsequent to t,o, the pupal population continues to increase 
until a maximum is reached at 7, due to increase in Zio. The writer is unable 
to say whether Z, is increasing or decreasing at this time. However, at pr, 
Equation 76 must hold. Immediately subsequent to p;, N's becomes negative 
by a continuation of the reasoning set forth in the above paragraph. 

During the period subsequent to ¢, the protective influence of the pupae 
and the progressive diminution in Z, produce a slight increase in the rate 
of increase of larvae in general. 
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To return to a consideration of .: the reversal of the sign of the derivative 
N’, is due to transformation of the predatory but non-egg-laying immature 


adults into predatory but egg-laying mature adults. 


Thus by the operation 
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Fic. 5. Smoothed curves for population growth of T. confusum, Duv. at 22°C. 


of Zi, the function descriptive of the rate of this transformation, N, is greatly 
increased, and thus ReN, is increased. This results in an increase in Ns, 
so that it becomes greater than zero, and the egg population rises to reach a 
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maximum again at py. 

This further arrest in the 
increase in the number of 
eggs is due to the with- 
drawal of pupae as an 
alternative food for the 
predators subsequent to pp. 
A similar reduction in the 
numbers of larvae is seen 
immediately subsequent to 


Subsequent to p», by the 
continued operation of this 
factor, N’: again becomes 
negative, and the egg popu- 
lation falls to another mini- 
mumat py, where, by reason 
of the continued increase 
in Ny, the sign is again 
reversed. 

Subsequent to pio, the 
egg production increases 
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enormously, and here again the protective influence of an alternative food is 
seen, as before any sensible increase in the rate of hatching can take place, an 
increase in the larval population is seen at py. 

This effect upon the lar- 


vae is in turn reflected in = bE Bea eA | | | | | {i 
a slight reduction in N’’: as eee Tb Pde g 


seen at Pi (p12 is not clear- 
ly defined at 22°C. owing 
to slight larval activity at 
this temperature) following 
which the egg population 
continues to increase. 
During this enormous 
rise in the egg population, 
although C’’s of the form 
C’,z are very great, not 


all eggs can be consumed las zy re 
a Pe 
by the adults and the com- 
paratively small larval bie ere 
populations before the eggs ILE LELFELLLL 


have reached an age I>. el 
- 2s ost 


Hence some hatching does Boe SESE et 


occur to cause an increase 

in the larval population at 

pu. This increase in turn “1G. 7. Smoothed curves for population growth of 
so augments the predatory T. confusum, Duv. at 27°C. 

population as to again 

reverse the sign of N’’; and decrease the egg population-subsequent to px. 

This periodic fluctuation in the numbers of eggs and larvae is believed to 
continue indefinitely. In the case of two populations grown at the University 
of Minnesota, owing to an infestation of intestinal parasites (Gregarina), pupae 
were unable to transform to immature adults, and slowly withered away a 
few days after pupation. In these two cases the egg and larval populations 
fluctuated regularly and repeatedly, with a period roughly equal to 2 r; 
while the adult population slowly dwindled away as explained on page 661. 
In such a case, the population as a whole would in time die out completely, 
unless the parasites decreased sufficiently to allow a few new adults to emerge. 
This special case will be described in a later paper. 

Meanwhile, the adult population has risen to a maximum where Z,)=0 and 
from then on, slowly declines, while Z,) remains practically zero, as explained 
on page 661.. It will be noticed that as Ny, decreases, the amplitude of the 
pupal fluctuations becomes greater. While Ni, remains greater than Y 
(see page 662) few or no pupae emerge, despite the brief increases in the pupal 
population, but as soon as N,, falls below Y, as explained above, a new pro- 
duction of adults begins the cycle again. 
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The writer does not believe that two identical cycles can ever follow one 
another in finite time, as the population at any instant is a function of its 
whole past history. It is highly probable however that succeeding cycles 
approach some fixed form, reaching the same after an infinite time. 
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